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Abstract 

We develop a general approach to nonequilibrium nanostructures formed by one-dimensional channels cou- 
pled by tunnel junctions and/or by impurity scattering. The formalism is based on nonequilibrium version 
of functional bosonization. A central role in this approach is played by the Keldysh action that has a form 
reminiscent of the theory of full counting statistics. To proceed with evaluation of physical observables, we 
assume the weak-tunneling regime and develop a real-time instanton method. A detailed exposition of the 
formalism is supplemented by two important applications: (i) tunneling into a biased Luttinger liquid with 
an impurity, and (ii) quantum-Hall Fabry-Perot interferometry. 
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1. Introduction 

Non-equilibrium electronic phenomena in nanostructures represent one of central directions of the modern 
condensed matter physics [l|, Q • Advances of nano fabrication have allowed researchers to explore experi- 
mentally transport properties of a great variety of nanodevices. Many remarkable phenomena have been 
observed in far-from-equilibrium regimes, i.e. for sufficiently large applied bias voltages. 

A particularly important class of nanostructures is represented by coupled one-dimensional (ID) channels 
(quantum wires). The coupling may be due to tunneling between the wires or due to impurity- induced 
backscattering. Realizations of ID elements that may serve as building blocks of such networks include, in 
particular, semiconducting and metallic quantum wires, carbon nanotubes, and quantum Hall edge states. 

The standard analytical approach to interacting ID systems (Luttinger liquids) is the bosonization Q. 
Recently, a nonequilibrium generalization of the bosonization framework |J| was developed for setups where 
a nonequilibrium fcrmionic distribution is created outside of the interacting region and "injected" into the 
Luttinger liquid. We will focus here on a more complicated situation when the tunneling or the impurity 
backscattering takes place inside the interacting part of the system. Such coupling terms represent in general 
a very serious complication for the full bosonization approach, and we are not aware of any way to solve the 
problem exactly. We choose instead an alternative route based on the functional bosonization formalism 
[5] that retains both fermionic and bosonic degrees of freedom. Combining the functional bosonization idea 
with the Keldysh nonequilibrium framework, we derive Keldysh action for the considered class of problems. 
This action has a structure reminiscent of that of the generating function of the full counting statistics [H, 0] • 
Our action generalizes that of Ref. [8[ where a local scatterer under nonequilibrium conditions was explored. 

To proceed with evaluation of the functional integral, we assume a weak tunneling and develop a real- 
time instanton (saddle-point) method. This allows us to determine Keldysh Green functions characterizing 
physical observables under interest (tunneling density of states, distribution functions, current-voltage char- 
acteristics, etc.). 

The goal of this article is twofold. First, we give a detailed exposition of the theoretical framework. 
Second, to illustrate the method, we present application of this approach to two important problems: (i) 
tunneling into a biased Luttinger liquid with an impurity, and (ii) quantum-Hall Fabry-Perot interferometry. 
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The structure of the paper is as follows. Our general formalism is presented in Sec. [5J In Sec. |3]we use 
the approach to calculate the tunneling density of states of a biased Luttinger liquid with an impurity. In 
Seed] we apply the method to explore a nonequilibrium quantum-Hall Fabry-Perot interferometer. Section 
[5] includes a summary of our results and a discussion of prospective research directions. Technical details 
are presented in five Appendices. 

Some of the results of this article have been previously published in a concise form in short communica- 
tions 



2. General Framework 



propagate along the channels, denoted by lower Greek index /i, with constant velocity 



2.1. Model and Functional Bosonization 

Let us consider a general model of the ballistic conductor, which can be represented as a network of 
one-dimensional (ID) chiral channels and point scatterers, as shown in Fig. [1] It is assumed that electrons 

from source to 

^ and x®, respectively. In the physical world such channels are 
realized by quantum Hall (QH) edge states or right/left-moving ID states in quantum wires. At point 
scatterers, denoted by Latin index j, instantaneous tunneling between different channels occurs, which is 
described by the scattering matrix s- 7 . Typical examples of scatterers are quantum point contacts (QPCs) 
or impurities in nanowires. Somewhat less trivial type of scatterer is a multi-terminal junction that can be 
realized by a quantum dot under assumption that its Thouless energy is well above all typical energy scales 
of the problem such as the temperature and the voltage. 

Albeit quite simple, our quantum-wire network model covers a broad class of mesoscopic ballistic de- 
vices, including QH interferometers and quantum wire junctions (See Fig. Q}. We note also that the im- 
portance of network models has been well appreciated in the context of the integer QH effect, where the 
Chalker-Coddington network model 11 serves as a highly useful starting point for numerical and analytical 
investigation of the QH transition. 




Figure 1: Two possible realizations of our model: (a) Mach-Zehnder quantum Hall interferometer and (b) junction of three 
quantum wires. Channels are represented by solid lines, point scatterers by white circles, reservoirs by boxes. Arrows indicate 
directions of motion. 



To warm up, we briefly recall the construction of the bosonized Keldysh action in the absence of tun- 
neling, i.e. when all scattering matrices are trivial, s J = 1, and thus different chiral channel are fully 
disconnected from each other. We also require that each chiral channel in the absence of tunneling is con- 
nected to one source and one drain reservoir (rather than forms a loop). Following the logic of the Keldysh 
formalism [HQ, all field degrees of freedom are doubled: upper Greek indices a = —/+ denote fields tp a 
on the forward/backward branch of Keldysh's time contour C; integration along C is to be understood as 
J c dt'A(t') = Jdt'(A-(t')-A+(t')). 

The fcrmionic action of the ID fermions reads 



A" 



dt dx^^idt + uvA)VV 



77 ^ / dtdxdx' U^ u (x,x')Q fi (x)g v (x'). 
.... J c 



Electron-electron interaction is taken into account by potential U^v, charge density in channel /x is g^ix) ~ 
if) ^{x + 0)rp ^(x) with Grassmannian fields ^^(x), ip^x). Spatial integration extends along the corresponding 
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channels: x^ < x < x® . To proceed we apply the method of functional bosonization [B[ and decouple 
interaction via Hubbard-Stratonovich (HS) transformation, introducing the bosonic field ip^: 

Ao[ip,-ip,(p} = / dtdxrf^idt + ivtidx- tp M )%l> ti + -^>2 / dtdxdx' ip M (x)U~J-{x > x')<p v (x r ). 

fj, J c [iv •' C 

The fact that we consider ID chiral channels enables us to decouple VvVV< ano - Vn by a gauge transfor- 
mation ip^ — > e* e ? requiring 

(9t + v li d x )9* = -<p* (1) 

While resolving this gauge condition, one should properly take the Keldysh structure into account, which 
yields 

with £ = (x,t). Here the blocks of the particle- hole propagator Z?o M satisfy the relations 

(d t + v^d x )D T / f (£,0 = ±6(t-z'), 

{dt + v^D^t,?) = (3) 

and therefore Eq. ^ indeed solves the gauge condition (fTJ) . In the frequency-momentum representation the 
bare retarded/advanced particle-hole propagator Dq^ in channel fi is given by 

D r Q /a (u;,q) = , W± = W ±iO. (4) 

The Keldysh propagator depends on the nonequilibrium state of the system. In what follows, we 
consider the zero temperature limit. Under this assumption, the electron distribution functions, = fx, 
f^ = 1 — fx, of source reservoirs A are completely characterized by the applied bias n\ — eV\, 

f}{t,t') = f}{t-t')=e^^f}{t-t') with = zpJL^-L-, (5) 

where fo{t) is the real-time representation of the Fermi distribution function and a is a short-time cutoff. 
The components of the bare particle-hole propagators are then 

D%[x,t) = Iv^n^t-x/vJ, D^ f (x,t)=0(±t)D> fl (x,t)+9( T t)D< fl (x,t), (6) 

with equilibrium Bose function rig(t) = —i/2n(t =p ia). Further, the Keldysh particle- hole propagator is 
given by the equilibrium relation 

q) = (£>S„(w, q) - D^{u, q)) sgn(w). (7) 

The gauge transformation has a non-trivial Jacobian contributing to the action. According to the 
Dzyaloshinskii-Larkin theorem [l2| , this contribution is Gaussian, and the new action reads 

m<p\ = ^E / *z*z'Mt)v- i &e)Me) - E / ^(omo , (») 

with excess charge go a = n\/ (2~k\v\\), effective interaction 

vrf-m, a = u- v \x, x')s(t - o - ^n,(e - r), 
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polarization operator n^(£) = — «Gq^(£)Gq^(— £), and free Green's function 

After the standard rotation in the Keldysh space 0, Q and the transformation into frequency-momentum 
representation, one obtains the retarded/advanced components of the ID polarization operator 

r I (10) 

With the action -4o[y] being Gaussian, the respective average value (w(C))o and ^ ne corre l a t° r of 
the fluctuations <fyv(£) = - (p(£))o are simply given by (^(0)q = E„ Jc d £' $')Q0v(?) and 

From @, or symbolically 6^ = —D^^ip^, one obtains for the correlator iDe )jUI/ (£, £') = (<5O AI (^)50 y (^')) 
the relation 

£>e, M „ = -D^V^D^. (11) 

After having reviewed the results for the action in the absence of tunneling, we are ready to consider the 
case of a connected quantum network. This will be the subject of Sec. [ 



2.2. Keldysh Action 

With all preliminaries we are now in a position to formulate the Keldysh action of the connected quantum 
network when at least one node is characterized by a non-trivial scattering matrix s J ^ 1. For the case of 
a single compact scatterer such an action has been constructed in Ref. [8j with the use of nonequilibrium 
Green's function method. The result bears connection with the solution of the problem of full counting 
statistics 0,0]. In our paper we generalize this approach to the situation with many scatterers. It turns out 
that the Keldysh action in this case can be written in terms of a full time-dependent single-particle scattering 
matrix (S-matrix) of the system in a given configuration of field ip a , which we denote S a — S[ip a ](t, t'), where 
a is the Keldysh index. Let us emphasize that the S'-matrix is non-local in time and takes different values 
on the forward and backward branches of the Keldysh contour. Our result reads: 



-ilnDet [l-f + S[<p+]*S[<p-]f]. (12) 

The last term in Eq. (|12[) is a functional determinant with respect to (real) time and channel indices. It is 
understood that / in the expression for the corresponding operator has the structure / M „(t, t') = S^f^t—t'), 
i.e., it is diagonal in channel representation, with f^(t) being the Fourier transform of the source distribution 
function connected to channel fi. The second term in (|12p represent anomalous contributions (related to the 
Schwinger anomaly) that have been already encountered before, see Eq. (J5J). They are most transparently 
written in the rotated Keldysh representation: ip c ^ = (tp~ ±ip+)/2, lfj a = (II J - ILj ± (II> - H<))/2 and 
integration is performed along the real time axis. 

A detailed derivation of the result <jT2j) , which employs ideas of Ref. 7], is presented in Appendix A. In 
view of the importance of this result, we give also its alternative proof (Appendix B), which follows closely 
the method of Ref. H. 

Construction of Scattering Matrix 

Let us now discuss how the S-matrix for the systems under consideration is constructed. The elements 
of S®(t,t') give the amplitude that a wave packet incident from source /i at time t' leaves the system at 

time t through the drain v. They are sums S"^(t, t') = A^ a (t, t') over all corresponding paths p formed 
by elements of the network. 

Figure [5] shows an exemplary path p through a network of channels and point scatterers. It consists of 
an alternating sequence of two types of processes: 
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Figure 2: Network built out of four channels fi, k, v, A and point scatterers 1, 2, 3, 4. Counting fields measure outgoing currents. 
Along interior parts of channels classical phases <^u/ K /„, e -g- due *° magnetic field, are accumulated. An exemplary path p is 
indicated by a dashed line. 



Electron propagation in the potential ip® between x l and x 3 , . leading to the accumulated phase 

0«(zV;t)(t) = -v- 1 f dx^(x,t-(aP-x)/v„). (13) 

J x z 

In addition to the Hubbard- Stratonovich field <p", there may be other time-depending phases ^ (e.g., 
induced by a magnetic field) contributing to Eq. ([TBI . Note that ^(x, x l ) satisfies the same differential 
equation (9 t + v iu ,d a; )'&"(x, x l ;t) — —(p*(x, t) as 9" (x, t), but has a simpler ("incomplete") Keldysh structure 
which involves only the retarded/advanced components of the bare particle- hole propagator Dq^. We refer 
to ^(xjX 1 ) as a kinematic phase. To take a finite flight time of electron between x % and x J into account, 
we introduce a "delay operator" 

/\^(x° ,x i ;t' ,t) = S(t' — t— (x° - x l )/Vfj_) (14) 
Then the amplitude of this process reads 

M^'^t,*') = e t§ ^ x3 ' xZ ^/\^x\x l -t,t'). (15) 
Indeed, consider the ID version of the Schrodinger equation on a directed link x l — > x 3 , 

idt^n = {-iv^dx + 0^)-0m- ( 16 ) 

Using the definition of the kinematic phase (|13p , this equation can be solved independently one each branch 
of the Keldysh contour yielding the relation 

M*3>t) = e"»^* i '' t) M*i>t-&-* i )M (17) 

which implies that the scattering matrix is given Eq. (|15[) . 

Scattering/tunneling at point scatterer j :. The amplitude of instantaneous scattering from channel \x to v is 

st ll {t',t) = si ll 8{t' -t). 

Passing the charge detector at drain fi, which is described by the counting field Xv As a special variant of 
a., our formalism includes the theory of full counting statistics. A counting field residing in the drain lead 
fj, measures the current flowing in that drain. The corresponding amplitude is 

X^t',t) = e- i ^S(t'-t). 
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Then the action (|12p enables us to express the cumulant generating function of the network as a functional 
integral over ip, 

Z(x)= fl['Dtf{x,t)exp{iA(<P,%}, (18) 
J n 

where vector \ combines counting fields Xp. m an drains. 

Finally, the amplitude A^ a of a path p is the path-ordered (real) time convolution ("latest to the left") 
of the amplitudes of its constituent processes. As an example, the amplitude indicated by the dashed line 
in Fig. breads 

A^(M') - [K M a v (x°,xl) s 2 VK M a K (xlxi) M«{xl,xD) (f,t') 

= 8(t-t'-r) e -^+^ slA^ exp{iC(^,^;«)+^«(^ 5 4;*-r 3 )+^(^,^;i-T3-r 2 )} , 

where 71/2/3 denote the flight times of the subpaths x^ — > x\, x\ — > x 2 K: x 1 , — > x® , and r = T\ + T2 + T3 is 
the total flight time. 

2.3. Weak Tunneling Expansion 

Due to the complex temporal behavior of the scattering matrix analytical evaluation of the functional 
determinant (|12|) is not feasible in general. An approximate treatment is possible if a weak tunneling at 
the point scatterers is assumed (i.e. the scattering matrix s^,„ close to S utIi ), and the ultimate goal of 
this section is the expansion of the action in the tunneling strength. Since in the absence of tunneling 
the network is described by the Gaussian action one can introduce the tunneling action At [tp], so 
that A[tp] = Ao[(f] + At [<p], where the expansion of At[(fi] starts from second-order terms with respect to the 
tunneling amplitudes at the point scatterers. In Appendix |Appcndix C we show that an exact representation 



of At [y] is given in terms of a modified ( "regularized" ) functional determinant 



At\ip] = —i In Dot 



i-f + s+^s-f 



(19) 



The new, "regularized" scattering matrix S here is constructed similarly to S. Each of its elements 

S^fat') — J2 P -^% is a sum over the same paths p which contribute to S^^{f,t') and connect the 

source \i with the drain v. Full and regularized amplitudes, A^ and A$[ respectively, differ in the partial 
amplitudes assigned to the elementary processes a. and b. which constitute a path p: 

Propagation between x l and x J .. Only the time delay is taken into account: 

M^,x i ) = A tl (x^x i ), 

while phase accumulation is shifted to 

Tunneling at point scatterers j.. The off-diagonal tunneling amplitudes become "dressed" by tunneling 
phases $S„(xM) = 6" (x^t) - 6°(4,i): 

^(i,0 = e < *^(- , .*V^(t-f). 

The phases O^ = —Do^ip^ are defined as in Sect. [2TT1 and can be modified by additional time-independent 
phase contributions due to e.g. magnetic or counting fields as follows. If the additional phase accumulated 
by an electron which propagates along a channel \i from a position x to the drain lead /1 is denoted as 
^.out( x )) then the phase 6^ is modified according to 

e«(x,t)^Q«(x,t)-^. ou t(x)- 
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In our previous example, Fig. [5J the regularized scattering amplitudes read 

Sll(t,t') = e^M^^-^WM-W-Hf^-x^ga S(t-t'), 



-lot 



(t,t 



5(t-t'). 



The regularized scattering matrix becomes trivial in the "clean" limit, S = 1, since all effects of interaction 
are now contained in the phases of the off-diagonal elements of the regularized scattering matrices §1^ of 
connectors. Thus Eq. (|19j) can be expanded in (even) powers of the tunneling amplitudes: 



AM 



iV-Tr 



i - s + ^s- ) / 



We are now going to elaborate on the second-order terms in this series. 
Second Order Expansion 

We introduce a notation A = S~ . Up to third order corrections in the tunneling amplitudes [that 
we denote as C(tun 3 )] the tunneling action is 



At W\ =iTr 



=iTr 



(t-A)f + \{t-A)f{t-A)f 



(1 A)^^f^ + 2 ^ ] Afiu fi/Ai/fiffi 



0(tun 3 ) 
+ 0(tun 3 ). 



(20) 



In the last expression, the trace is only taken with respect to time. To reduce the tunneling action to this 
form, we used (1 — -4.) w = 0(tun 2 ) and A V{1 = O(tun). 

It can be shown (see Appendix |Appcndix D| for a detailed derivation) that At acquires contributions 
from paths which start in a certain source reservoir, evolve forward and backward in time, undergoing in 
total exactly 2 tunneling events, and eventually returning to the original source. Such paths involve exactly 
2 different channels, /i and v. Thus we can classify all paths according to the pair (p, v) of channels /i =/= v 
and the pair of scatterers (possibly equal) at which the tunneling takes place: /i — > v at i and v — ¥ [i at 
j. Of course, the class (ij; y,v) coincides with the class v\i). The second order expansion of the tunneling 
action then is a sum over these classes: 



AM = -i y d£ i di 2 (e-**^.**) 



-**£,0«:Vl) 




<i2 — ti — tl-, where the tunneling polarization operators are given by 



ngJ.W = [e(±t)u>.^(t) + e(Tt)n<.^(t)] xm 



T ■ — T - 
1 fi.m ' fj,.m 



T t:u\ 



(21) 

(22) 
(23) 



where r^ in is the flight time from the source to the scatterer k along a channel A, r^ in = (x\ — xf)/v\. 

We have also taken into account counting fields in the drain leads (which are not contained in n^/J^) and 
classical phases, 

A#f„ = ^.out(^) - </V.out(£ l ) - <pv.aut(x j ) + AI .out(^)- 



In the case i = j we will also use the convention 



n^W + n<^(i)] 



(24) 



The comparison of this expression with the Eq. (|23[) shows that they differ from each other by the singular 
term proportional to sign(t)h(t)5(t) — ir5 2 (t), where we put hit) — f^{t) — }^{t). It gives some constant 

(albeit infinite) contribution to the tunneling action (|2T|) and therefore both representation for TlJJ^ are 
equivalent. 

2.4- Real-Time Instanton Method 

On the level of the second order expansion, the action At [<p] is expressed in terms of the tunneling phases 
3>, which are linear functionals of cp: 

*M"(0 = Y, f tfVw\{S,?)<P\(]?), ^;A = VM-W»A. (25) 
A J C 

The action At[<p] is non-Gaussian in $ and, in fact, is quite similar to the Ambegaokar-Eckern-Schon (AES) 
action (13J. The difference is that the kernel appearing in Eq. (|2ip is not only non-local in time (as in the 
case of AES) but in general is non-local in space as well. In view of the non-Gaussian character of the action 
an exact evaluation of physical quantities does not seem feasible in general. For this reason, we will use a 
saddle-point approximation that catches correctly the relevant interaction-induced physics, including both 
the renormalization and the dephasing phenomena. 

To explain the idea of the method, let us consider some physical quantity 0[ip] = Oae lAj ^ v \ where 
»4j[v] = — X)u Je ^£ JfiiQVuiO i s a linear functional of tp, and the prefactor Oq is independent on p. 
Important examples, which are treated in the next sections, include the electronic Green's function and the 
current. The quantum average value of O is given by the functional integral 

(£%]) = / ^e^obl+a.M+a.M^^ ^ 



which we estimate in the semiclassical approximation [14]. In this case the path integral is contributed by 
the saddle-point trajectory ip** of the full action .4o[<p] + .At[<£] +_4j[</j] and quantum fluctuations around it. 
Here free and tunneling actions -Ao[^L A*M are given by Eqs. © and (|2ip . respectively. In the limit of weak 
tunneling between chiral channels the saddle-point trajectory ( "instanton" ) can be found approximately by 
requiring that it minimizes the Gaussian contributions to the action Aq [<p] + Aj [p] , which gives 

MO = E / d £' + Mi'))- (27) 

As shown in Appendix |Appcndix E| under such an approximation Eq. (|26[) simplifies to 

(0[cp]) = (e^M} e'&teJOo, (28) 



with 



{j Aj[v] ) Q = CX P {i (M<P\)o \ [((A/M - (^M) ) 2 )J ] , (29) 



v [d tl dt 2 r e -»i-c 4 A) e-^U'**)) ( n Jr ( !; 2 \ '^n'f 



A t [ip*] = —i 



(30) 



where t%2 = t\ — <2 and (. . .}q denotes averaging with respect to >4oM- We have introduced the instanton 
phases = V^tp* and the renormalized tunneling polarization operators 

fi$U(*i - h) = j{nilM^-M-t*)-D^ m ) n ^ (ti _ hl (31) 
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obtained by dressing of the bare tunneling polarization operators by phase-phase correlators 

The meaning of Eq. pip is that quantum fluctuations of tunneling phases renormalize the temporal depen- 
dence of tunneling polarization operators which lead to non-trivial (usually power-law) energy-dependence 
of tunneling coefficients. 

In the next two sections we consider two important applications of our general approach. 



3. Tunneling Density of States of Luttinger Liquid with Single Impurity 



In this section we show how our formalism can be applied to the evaluation of the tunneling density of 
states (TDOS) of a nonequilibrium quantum wire containing a single impurity [9j], as depicted in Fig- [31 On 
the experimental side, the interest to such a theoretical study is motivated by the recent development of the 
nonequilibrium tunneling spectroscopy of ID nanostructures, including carbon nanotubes [l5| and quantum 
Hall edges [H, Hill. 



3.1. Model and Results 




Figure 3: Tunneling experiment with a voltage-biased quantum wire. 

The wire is modeled as a network of two channels: right- (left-) moving electrons, denoted with indices 
fi = +(— ). Backscattering, i.e. tunneling between the two channels, occurs at the impurity (at x — 0) which 
is considered as a point scatterer with scattering matrix s. Non-equilibrium conditions are established by 
biasing the right reservoir with respect to the left one: eV+ = eV > 0, eV- = 0. 

The interplay of interaction, nonequilibrium, and impurity scattering can be studied by placing a con- 
ducting tip held at a voltage Vti p near a position x of the wire (we can assume without loss of generality that 
x > 0) and measuring the current /tun between between tip and wire. Let us assume that the tip can be 
described in terms of fermionic quasiparticles with density of states ^tip(e) and distribution function /ti P (e), 
as is the case, e.g., in the absence of interaction within the tip. If the coupling |£| 2 between the tip and the 
wire is weak, a simple perturbative expansion yields the tunneling current 




The "rates" for tunneling into/out of the /i-channel of the wire are defined as 




Measurement of a differential tunneling conductance in the limit of a small temperature gives the access to 
the tunneling density of states, since 

« r>{e) + T<(e)=^(e) = --ImGl(x,x,e). (33) 

Vt ip =e 7F 



dV 
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In the absence of interaction, the rates would simplify to rjj(e) = //f(e)^ M (e) with the distribution functions 
/< (e) = / M (e), fu(e) = 1 — /^(e) and the density of states v^{e) in the channel \i. Evaluation of tunneling 
rates and of the TDOS of the interacting problem is the goal of this section. 

We consider a spinless LL model with a point-like repulsive interaction, t/ M „(x,x') = VoS(x — x'), that 
does not discriminate between different channels. The interaction strength in the LL model is conventionally 
characterized by the constant K = (1 + Vq/ttvf)^ 1 ^ 2 ■ The free electron spectrum is linearized around the 
Fermi points, which requires the introduction of a high-energy cutoff A ~ Ep (which is of the order of the 
bandwidth). In the absence of backscattering, the tunneling rates exhibit the well-known zero bias anomaly, 
i.e. a power-law suppression near the Fermi edges, 



eVu 



A 



2l 



rj(e) = -r(i + 27)- 1 x 0(±( C - evy) 

where Vq = (2ttvp)^ 1 is the non-interacting density of states, and the exponent 7 is given by 

{l-Kf 
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4K 



(34) 



(35) 



As is shown in the remainder of this section the tunneling rates change considerably upon including the 
impurity. For eV > the rates are given by 



r|(e) = ±^ (jp) 27 r(-2 7 ) Im [( T z^ + WeF)(±l) 2 ^(-2 7 , 1 - 2 7 + 2S M , -1 - zj] 



where we have introduced the following notations: 
T{2K) 



C± 



, 5+ = {l- K)/2, 5-=l/2-K, ZfJ 



T(l/2±K/2y 

We have also introduced the renormalized reflection coefficient 



(e-eV M + -r- 1 )/eV. 



R*{eV) 



T{2K) 



A 



eV 



2(1-K) 



and the nonequilibrium dephasing rate 



R*{eV) 



2 sin ir8+ 



eV\. 



(36) 



(37) 



(38) 



(39) 



The energy dependence of rates Tig (e) is shown in Fig.H) The main feature of these plots is that the tunneling 
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Figure 4: Tunneling rates for right- and left-movers with different interaction strengths: K = 0.3 (1), K = 0.5 (2), and K = 0.75 
(3). Fermi edge for right- (left-) movers is at eV+ = eV (eV— = 0). 



rates have split power-law singularities which are characterized by different exponents and are smeared by the 

10 



nonequilibrium dephasing rate 1/r^,. The main edges are located at the corresponding chemical potentials, 
i.e., at e = eV± in the case of right-/left-moving states, respectively, and are characterized by the exponent 
equal (in the considered weak-back-scattering regime) to its equilibrium value 7. The formation of the 
second (side) edge due to scattering off the impurity occurs at e = e(V± — V). If the interaction is repulsive 
(K < 1) then the corresponding exponent 2(7 — <5_) for left-moving electrons is always positive, hence the 
correction at the side edge e = is smooth. For right-movers in the case of not too strong interaction, 
K > 1/3, the nonequilibrium exponent 2(7 — <5 + ) is negative, yielding a resonance in tunneling at the side 
edge e = 0. 

The presence of side edges in the tunneling rates can be understood in the following way. Inelastic 
electron backscattering at the impurity at point x — induces the emission of real nonequilibrium plasmons 
with typical frequencies fku < eV, which in the non-dissipative LL can propagate to the distant point of 
tunneling x. As the result, inelastic tunneling with absorption or stimulated emission of these real plasmons 
become possible. For example, an electron tunneling into the right/left moving state of the LL with the 
energy e < eV± can accommodate itself above the corresponding Fermi energy (eV±) by picking up the 
quantum hu > |e| from the nonequilibrium plasmon bath. Since the energy of out-of-equilibrium plasmons 
is limited by the applied voltage, one has a threshold: e > e(V± — V), which is developed into the power- law 
singularity typical for the LL. The singularity at the side edge of the tunneling rate-out describes the inverse 
processes: the inelastic tunneling from the LL accompanied by the stimulated emission of nonequilibrium 
plasmons with typical energy Huj ~ eV. Such side edge is pronounced in the case of right-moving states 
only, r$(e), and is not seen for the left-moving states, rf (e), since the associated exponent 2(7 — #_) is 
always positive in the latter case. 

Having announced the main results, we now turn to details of their derivation. 

3.2. Calculations 
3.2.1. Action 

Since interaction does not discriminate between channels [i,v — ±, it can be decoupled by a Hubbard- 
Stratonovich transformation introducing a single field if, i.e. (p+ — (p- — (p. For weak backscattering, i.e. 
weak tunneling between right- and left-moving states at the impurity, the action A[(p] — Ao[(p] + At[tp] is 
obtained according to Sect. 12.31 The free action © is 

m<p) = \ I dedevo^cWMO-E / d ^(0(eo+ + e -) 

J C ,. J c 

with £ = (x,t), non-local effective interaction V^ 1 ^— £') = V ~ 1 5(£ l — £') — n(£— £') and the total polarization 
operator II(£) = n + (£) + II_(£). Using (fT0|) one obtains for the retarded/advanced components of effective 
interaction 

V^,p) = V U}2 - VF X (40) 

with a plasmon velocity u = vf/K. (For a repulsive interaction K < 1, so that u > vp.) 

Since we are dealing with a single scatterer, the tunneling (or, equivalently, backscattering) action At, 
as given by (|2"Tj) . consists of one term [corresponding to class (11; H — )]: 

where § T (t) = 05(0, t) — 0ljl(O,i) is the tunneling phase evaluated at the impurity, x = 0. The phases 
are related to the Hubbard-Stratonovich field (p according to ([2]). The tunneling polarization operator 
11^ = IIii;_| is given by Eqs. (|2"2"|) . (|23l) . The components of the polarization operator read 

H$_(t) = -|r„| 2 /|(i)/ g (-t) = -M 2 e-^ \f^(t)} 2 , (42) 



n T + /f (t) = l[n>_(t) + n<_(t)]- (43) 
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where |ro| 2 = \sh | 2 = |sl + | 2 is the bare reflection coefficient. We note that in Ref. @ different notations 

were used: Ab instead of *4o[<^], A[ mp instead of -4t[y], and $ instead of 0. We also note that Eq. (|4"5)l is 
slightly different from Eq. (051) : however, this difference leads only to an additional constant contribution to 
the action that is physically irrelevant. 

3.2.2. Green's Functions in Instanton Approximation 

In order to find the tunneling rates, we represent the electron Green's function at the point of tunneling 
x > as a path integral over the field tp, 




Here, G fl (x,x;t,[(p}) denotes the Green's function for a given configuration of ip. It satisfies the Dyson 
equation with the spatially local self-energy 

S M M(x, a; ^ ) 0=^(^)5(xO(|ro|V/2)e i ' 1 *« 5 _ M (t-Oe-^ (t ' ) , 

where are the quasiclassical Green's functions of the source reservoirs. Solving the Dyson equation to 
the first order in |r | 2 , we get 

Gf(x,x;t)=g^+g^, (44) 

where 

gap = ^m G f^ £; t) e -<(^ j 

<^E^fMt 2 (45) 

7<5=T 

(e'^^G^i.O;^ tJeWMglfa ~ i 2 ) e -^^G^(0, x; h - t) e - <e 2< s '°>) . (46) 
Here Gq^ are the Green's functions of free electrons; in particular, 

G^(x,t) = ±/*(t - »x/v F )/iv F . (47) 

All averages (...) in Eq. (|4"6")l are taken with respect to the action -4o[<^] + -4t[y>]- They are of the form 
(|26[) and can be evaluated with the real-time instanton method described in Sect. 12.41 In this approximation 
the first term in Eq. (|44|) reads 

S5S « e iMv * ] x (e m i^G%{x,x-t)e~ l ^^) a . (48) 
The second factor here is the full Green's function of a clean LL, 

a> - t\ (s,t)-ef (x,o)] 2 ) „> ,_ _ ^ a 27 * 1 

G^(x,x;t) = e "Um( )J ;„ G <^ X]T) = ± __ e ^ (g ± . (49) 

The first factor in Eq. (|48p gives dephasing corrections due to the interplay of tunneling and interaction. The 
instanton action ^[t^,] is defined in (|30p and obtained by substituting the dressed polarization operators 
into (14"T1) . The instanton phase ip* is generated by the source i^4j[y] = iQ®(x,t) — i@P(x,0), 

^(t) = (^(t)) -D^(t-i,-x)+D^(t,-x) with D&fi,x)=Dl*_^t,x)-Dl* +ll (t,x), (50) 

and depends on Keldysh indices a, (3, the direction /x of the tunneling electron, as well as the time and the 
position of tunneling t, x. Since the average value (3> =F (i)) does not depend on the Keldysh index qp and 
the time t it will drop out when the instanton phase is substituted into (14"T1) . Therefore it will be omitted 
in what follows. 
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3.2.3. Phase-phase Correlation Functions 

To lay the groundwork for all further calculations we compute the correlation functions of the phases 
M and $ = 0_ — 8+. With the bare particle-hole propagator the effective interaction (|4"0")l . and the 
relation (TTT|) . the retarded/advanced components of the correlators D®^ v can be easily evaluated in the 
(q, ui) representation: 



D 



r / a 



d 



r / a 



( w ,<7) 



2ir uj + [ivpq 

U) U) 

2tt 1 - K 2 r 



-/' 



1 + K 



1 



1 -K 



1 



1 



1 



u) AK 



AK q — /j,uj±/u 4K q + nuj±/u 2q — /jij±/vf 
1 1 



q — oj±/u q + uj±/u 



Transforming the above relations into the mixed space-frequency representation, we obtain 



D 



r/a 



2ni 



with 



c,;„ e 

= (1 ± K) 2 /(AK), c± 



9(t^x) c e 



— —ifiujx/u I 
ftp 



= (l-K 2 )/(4K). 

The Keldysh component of the phase correlator is given at zero temperature by 

D%^,x)=B{u){Dl^ v {w,x)-D a e ^ v {u,x)) = S ignu, (^(m) - D^(oj,x)) . 
Performing the Keldysh rotation, we then arrive at 

^IL>' x ) = \ ( D e,^(", x) ± (jr eiMV ( U , x) - Z>^>, x))) = ±0(±u>) (D r e ^(u>, x) - £>|, M >, x)) , 



(51) 
(52) 

(53) 
(54) 
(55) 



(56) 

In the real-time representation (x, t) these phase-phase correlation functions can be decomposed into the 
plasmon (moving with velocity u) and particle- hole (having velocity vp) contributions 



a0 



(57) 



where for a given velocity v the functions (t, x) read 



C$(t,x)= In 



i =F ia — x/f 



, C T J f {t,x)=\n 



=Fia signx/u 
t =p ia signx/u — x/v 



This follows from the Eqs. f|56[) after the Fourier transformation from u> to t with taking into account the 
high-energy cut-off. The functions (i, x) satisfy 

d t £^(t,x) = —(t^fia — x/v)" 1 , d t C^ T (t, x) = —(t^fia signx/u — x/u) _1 . 

It is worth mentioning that the appearance of both plasmon and particle-hole "light-cone" singularities in 
the phase-phase correlation function is a special feature of the functional bosonization approach. 
For the — $ phase-phase correlators we then obtain: 



iD%Ut,x) ee (5^(t,x)6QP (0,0)) = i \D$Jt,x) - D$ + Jt,x) 



p£f u (t,x)~q£^ u (t,x)~£f VF (t,x) 



The correlation function of the tunneling phases $(i) = &(t,0) at the position of the impurity reads 



iDf(t) ee (5$ a (t)5$P(0)) = i lim \Df e _(t, x) - Df @+ (t,x)] = -2(1 - K)C% p {t) 



where 



4W = m^, 



4 /f (*) = J 



hi- 



t — ia 



In- 



t + ia 



(58) 
(59) 

(60) 
(61) 
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3.2.4- Instanton Action 

The correlators obtained in the previous section reduce the "dressed" tunneling polarization opera- 
tors (|3"Tj) to the form 



n$_(*) = -M 



-ieVt 



(2na) 2 



a ± it 



2 A' 



(62) 



or in the frequency representation 



R*(eV) 
2ir 1 



(±(w - eV)) 



2if-l 



|eF| 



where we used the definition for the renormalized reflection coefficient i?* . With the mixed phase-phase 
correlation function (|58[) at hand we are also in a position to write down the instanton trajectories (|50[) , 



(t) = /x {phi 



t ±ia + \ixju 



t ± ia — t + /j,x / u 
a — ifi{t + fix/u) 



q In 



a — i/3(t — jjix/u) 



a — ia(t — t + fix/u) 



gin 



ia(£ — i — fix/u) 
t±ia—[ix/u 



t ± ia — t — fix/u 



- In 
-In 



t ± ia + fjLx/vF 
t ± ia — i + fix/vF 
a — i/3(t + /xx/vf) 



a — ia(t — t + (ix/vp) 



[ix > 0, 
, [ix < 0, 



where we have introduced p = (1 + -KT)/2 and g = (1— K)/2. It is worth emphasizing that these instantons 
represent non-classical solutions in the sense of the Keldysh nonequilibrium theory: the phases are 
different on the upper and lower time contours, so that the quantum part is non-zero, ^ 0. Because 

of this the corresponding tunneling action y^^*], which we are going to evaluate, is non-zero. 

To exemplify the evaluation of the instanton action _4.{[<p*], let us consider the case of tunneling into/out 
of a right-moving state with the tip being placed on the right from the impurity, u — + and x > 0. The 
phase factor is 

e^ (t) =^(f) K !(f)4(t) 



with 
/£(<) 



t ± ia + x/u 



v t Ha — t + x j u 
The instanton action reads 



a — ia(t — t — x/u) 
a — i/3(t — x/u) 



«?(*) 



t ± ia — t + x/vp 
t ± ia + x/vf 



(63) 



(64) 



(65) 



%M<P*\ = & kh&uft;>_(h-U) K%>(t at U) K. c J(t 3 ,U) KtffaU) 

with K£"(t 3 ,U) = i4{UT 1 kI{U) for (7 = +,-,0. 

Since the polarization factor n_| (£3 — ^4) comes with the factor e - ieVr (*3-*4) ^ e integral is dominated by the 

region \t$ — £4 ] < |ey| _1 . Furthermore, important contributions are expected to come from regions around 
the singularities of the phase factors, i.e. £3,1^4 ~ —x/u, i — x/u for /C+, £3,44 ~ x/u, i+ x/u for /C_, and 
t3,ti ~ —x/vf, t — x/vp for /Co- These regions in the £3 — i4-plane are sketched in Fig. [5] 
We will assume that the singularities are well separated, which imposes the condition 



\t\,\eV\- l ^(l-K)\x\/v F) 



(66) 



For interaction strength of order unity, the dephasing time t v (which governs the relevant t) is ~ |ey| _1 , 
so that the conditions (|66p reduce simply to leV^" 1 -C \x\/vp. This condition, implying a sufficiently large 
voltage and/or tip-to-impurity distance, can be easily satisfied. 

Far from the singularities, the phase factors become trivial, K. a (tz,ti) — > 1, so that the integral (|65l) 
approximately splits into 

iAM t»X[K+]+X[K-]+I[lCo], (67) 
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-x/v F -x/u 




~(eV)- 




x/u I 



Figure 5: Regions in the t$ — t4-plane providing dominant contributions to the integral for the instanton action, Eq. 113,2. il l 



with 



x[K a \ = J2 to /dt 3 dt 4 n c + I '_(t 3 - h) (K ( J>(t 3 ,u) - 1) . 



(68) 



C,V=T 



We have added —1 to the phase factor YL G in Eq. (1551) to make the convergence manifest. This does not 
change the value of the integral in view 

of E C r,to&+-{t) =0. For the very same reason, the independence 



of K_ (£3,^4) of the Keldysh indices £, 77 implies I[/C_] = 0. 

According to Eqs. ((MJ), ([5X5]) . /C+, /C have the form JC ir >(t 3 ,U) = k c (ts)' 1 n n (t 4 ) with kt^) = 
( t-ttSa ) an d some exponent r > 0. Therefore, X[/C] is dominated by the regions \ta — < |e V"| 1 

(singularity of 11^ ) and — t a |,|f 4 — tb\ < |eV] 1 (singularity of K.) which determine the long-time 

asymptotics \eVi\ 1. This yields 



I[K] « fdtJ2toH+-{t) fdT (K C (r)- 1 /fi"(T) - 1 



i?»(eT/) 



|eV(t 6 



2tt 







—2ixir sign [eV(t b — t a )] 



r(2K) 



ieV(t b -t a ) 



ieV{t b - to)] 



2(r-K) 



(69) 



Assuming for definiteness x > and considering first right-movers, fi = 
action L4 t [</>*] — I[/C+] +I[/Cq]- The first contribution here, 



we thus obtain for the instanton 



l[K + ] = -f-iR*(eV) S ^eVi 

2T<p lit 



C+R*(eV)e leVt [ieVt\ 2q 1 C+ = T{2K)/T{pf 



(70) 



encodes effects of real plasmons on tunneling which are generated because of backscattering off the impurity. 
One of such effects is the shot-noise, which is represented by the first term in Eq. ([70]) . This term is negative 
and linear in time t and thus accounts for dephasing with the rate (|59")l . The second term in Eq. ([TDl represents 
a perturbatively small renormalization of bias voltage and will be neglected in the following. The third term 
in Eq. (jTO")) is subleading as compared to the first one and we will treat it perturbatively in i?* . It shows an 
oscillatory behavior accounting for an energy transfer ~ eV between the nonequilibrium bath of plasmons 
and a tunneling electron. We will return to this point when discussing the tunneling rates. 

We turn now to the contribution Z[/Cq]. In this case r = 1, so that the first term in (1691) vanishes. Thus, 
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the the electron-hole pair contribution reads 

I[/C ] = E^n c + "_(f) [dtse-^nifo)- 1 /di 4 e ieyt ^(i 4 ) (71) 

*—* eV=0 J J 

Oi 

= R*(eV)T(2K)e- leVi [-ieVt\ 2(1 ~ K) . (72) 

It is oscillatory and can be again treated perturbatively. In contrast to the plasmon contribution, however, 
it seemingly corresponds to an energy transfer of ~ — eV and does not have a clear physical interpretation. 
We will see below that this term is an artifact of functional bosonization which will be canceled by the Born 
correction Q"j^ . In total, we have 



e iA t [<p,] w e -\i\/2r v fi + C+ R,( e V) e ieVi [ieVt\ 2q +X[K ]) for x > 0, fi = +. (73) 
The same considerations can be applied to left-movers, with the result 

e iA t [v«] ~ e -|*l/2T„ A + C-R*(eV)e ieVl [ieVt\ x ~ 2K \ , C_ = T{2K)/T{q) 2 , for x > 0, /i = -. (74) 
3.2.5. Bom Correction 

We evaluate the Born correction Q"^ to the Green's function, (14"4"]| . in leading order in |ro| 2 , which 
amounts to taking averages with respect to the clean action only. Then Wick's theorem yields 

£ l5 UdtzG^i^O^-t^ih-t^G^ (75) 

7<5=T 

with 

J^ 5 (x,t;h,h) = ^ e ^^(s.t) e iM^(t 1 ) e -^* s (t2) e -ie^(x,o)^ 

= e -i<[eS(s.*)-eS(e.o)] a ) o -i<[»''(ti)-»'(* a )] a > o ^(^(tx)-*^)) (76) 



and the instanton (|50[) . The appearance of the instanton makes the integral (|75[) quite similar to the 
instanton action and we will use an analogous approximations to deal with the time integrals. 

We have already seen that the instanton phase factor factorizes into three contributions (|63p — two gov- 
erned by plasmons and one by electron-hole pairs — and one might expect the integral (|75[) to split into three 
contributions in a way akin to the instanton action. However, it will turn out that the presence of the bare 
Green's functions Gq^(x, 0; i—ti) and Gq^(0, x;h — t) suppresses the plasmon contributions. Focusing again 
on /i = +, x > 0, we show that the remaining electron-hole pair term cancels I[/Cq] in (|73[) . 

The t\ , i^-dependent contributions to (|75p are 

G$(M;*-*i)s:i tf (ti-*2)G&(o^ 

= e- ieVl e ieV ^-^ x fj(h - t 2 ) x 4{t 2 )~ l K 5 + {t 2 )~ l k 5 _ (h)' 1 x 4(tiW + (tiW_(h). 
We combined terms with similar pole structure, defining 

gL(h~h) = e -\ L J /o^ (tx -h) 



eV=0 



4(h)- 1 = G° p + (-x,t 2 ) 4(h)- 1 , 

eV=0 

4(h) = G^(x,i-h) 4(h). 

eV—O 

All voltage dependence has been singled out in the phase factors explicitly. 

16 



Let us examine the pole structure: g1_ S (ti —£2) is divergent for t\ f=a £2- This is reminiscent of (£3— t^) 

in (1651) which preferred i 3 « i 4 . The plasmon contributions k + and k_ have been studied in the previous 
section where we already noted K±(t) — > 1 for t far from their singularities. This is no longer true for kq- 
Indeed, leaving the Keldysh indices and the corresponding short-time regularizations aside for a moment, 
we have 



G a+ (x,t) 



1 



1 



eV=0 



2-kvf t — x/vf 



and Ko(t) 



t — t + x/vp 
t + x/vf 



and hence, 



4(h)- 1 



2ttvf t% —t-\- x/vp 



1 



2nvp t\ + x/vf 



Similarly to the original phase factors Kq^) _1 , k o(^i)i these new ones have the poles at t\ ~ —x/vf, 
i 2 ~ i — x/vf', however, at variance with Kq^) -1 , «o(ti), the factors Kq^) -1 , ^o(^l) vanish far from the 
poles instead of converging to unity. Therefore, the poles of Kofa) -1 , &o(ii) are dominating the integral 
([75]). while the plasmonic poles give subleading contributions (suppressed by the factor Vpt/{1 — K)x <C 1). 
With only leading terms taken into account the integrals simplify to 



a/3 



2 



J2 ^e-^^-'^ [Gr + (xJ+x/vF)gl\-t)Gli(-x,t-x/vF) 

7,<5=T 



eV=0 



dt 2 e- ieVt * 4(h)- 1 dt ie teVt i K^ti) (77) 



For large times \t\ ^ \eV\ , the short-time regularization and thus the distinction between different Keldysh 



components becomes immaterial, e.g. IfL*_ 

Gq+(x, t + x/vF)gl S (—t)G 5 +(—x,i— x/v F ) 



r \ 2 e- ieVt /(2Tri) 2 , G^(0,t) = -l/(2-rTv F t), and one obtains 



\r \ 2 v F 



-ieVt 



eV=0 



1 1 



2itvf t 



i 1 



-ieVi 



= -G$(S,S;t) 



57 



e v=o 



Substituting this into (j77|) . taking into account the dressing of the Green functions (|49|) and polarization 
operators by phase factors and comparing with Eq. (1711) , we find 



^7 

-Go+(»,x;t) I[/C ] 



(78) 



The very same analysis can be performed for fj, = — . In this case, however, Kq does not depend on the 
Keldysh index 7. Because of X) 7 5 7^ 9+(ti ~ ^2) = the contribution is negligible. 
Concluding, we obtain 



0£f « e -l f l/ 2r ^ G^(x,x;t) (l + C + i?,(ey)e ieyf [iel/tf 9 +I[/C ]) , 
^f«e-l^ 2r " G£f (l + C_J2.(eV) e ie 



1-2 if 



G%l(x,x;t) 1[Kq] 



(79) 
(80) 
(81) 



In leading order in |r | 2 , i.e. neglecting dephasing corrections, cancels the Z[/Co] term of Gq+, as was 
stated in the end of Sec. 13.2.41 
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3.2.6. Tunneling Rates 

We are now ready to evaluate the components of the Keldysh Green function and, in particular, the 
tunneling density of states controlling the tunneling current between the system and the tip. Let us note 
that because of the tunneling actions (l?3"f and ([74")l the Green functions contain power-law terms [iel/i] 25 , 
r > 0, which have apparent branchcut singularities near t = 0. However, results ((73)) and l|74|) are only valid 
in the long-time limit \eVi\ 3> 1. A regularization which takes this into account and does not violate the 
> 1 * > 

symmetry relation iG~ji (x, x;t) — iG^ (x, x; — t)) is [1 + ieVt\ 2S , which yields the Green's functions 

Gf(x,x;t) = Gf{x,x-,t)e-W/ 2T * (l + C v ,R„{eV)e ieVl [\ + ieVt] 25 ^ (82) 
with C+ = T{2K)/T{p) 2 , C- =T(2K)/T{q) 2 , 5+ = (l-K)/2, 6^ = 1/2- K. (83) 

The tunneling rates are obtained by Fourier transformation to energy representation. Using the aforemen- 
tioned symmetry property of the Green's function, we obtain 



r«£(e) = ±— / dt (e iet iG^(x,x;t) + e-' et iG*(x,x; -t) ) = T -Im / dt e ltt G^{x,x;t) 
2tt J V / n J Q 



= ± — Im 



J<(e; eV,) + C»R*{eV) J s <(e; eV, - eV) (84) 



Jf(e;U) = ±ia^ / dte^- u+ ^ 2T ^ {1 + ieV^ 25 {a± it)~^ +1 \ 



with 



Since J^{e\ U) is an analytic function of all parameters e, <5, 7, £/, eV^ € C as long as Im (e — U+i/2T v ) > and 
ImeV < 0, we can consider here Re 6 > — 1, Re (27 + 1) < 1, ImeV < and Re [{e — U +i/2T lp )/eV] < and 
deduce all relevant cases by analytic continuation. Under these constraints, one can evaluate the integral by 
rotating the integration contour, t = —is/eV, < s < 00, into the complex plane, and put a — >• 0. Writing 
z = (e — U + i/2r v )/eV we get 

, 2 7 



poo 

J*(e; U) = ±— / dse zs (1 + *) M (±«/eV)- (27+1) 



eV\ 21 

±^-J r(-2 7 )*(-27, 1 - 2 7 + 26, -z) (85) 

with the confluent hypergeometric function \&(a, b, z). In the case 6 = one has ^(—27, 1 — 2j,z) = z 21 
and thus the term J§ in Eq. (|84l) yields the equilibrium zero-bias anomaly. Taking into account the 2nd 
(impurity) contribution we arrive at the results that have been presented in Sec. 13.11 [ see Eq. (f3l)|) ]. We 
note here that the result for the physical (real) voltage eV is obtained from Eq. ([551) by substituting there 
eV - iO. One can also note that at 5 ^ one gets *(-27, 1 - 27 + 26, z) — z 2(7_<5) at z 0. Therefore 
Eq. I|85p gives the impurity correction Ar M to tunneling rates which is singular at e = eV± — eV. Explicitly 
one has 

-<M~_p ^L_.v, ^.t/I 2 ^-^)^/ 27n5 M- £>eU±-eU 



AT<( 6 )«-^(eV)| e -eV ± + eVr™ sin j ^ Y< e V±- e V (86) 



in case of tunneling from the wire into the tip and 



AT>(e) ozR*(eV) 6{e - eV ± + eV)\e - eV± + eV\ 2{l S " ) (87) 



in case of tunneling from the tip into the wire. 
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4. Quantum Hall Fabry-Perot Interferometer 



In this section we study the role of Coulomb interaction in an electronic Fabry-Perot interferometer 
realized with chiral edge states in the integer QHE regime. Electronic Fabry-Perot (FPI) [26[ 28, 27, 2^] 
and Mach-Zchnder (MZI) [13, EH H3, [H, [Hi, [BH [H, [H3, [H, [H, S3, EU interferometers are analogues of 
the optical interferometers, where the chiral edge states play the role of light beams while quantum point 
contacts (QPCs) act as beam splitters. Electron interferometry provides a powerful tool for studying the 
quantum interference and dephasing in mesoscopic semiconductor devices. Another motivation behind these 
experimental efforts stems from the recent interest in topological quantum computations, which propose to 
exploit the non-Abelian anyons in the fractional QHE regime [l9j ]. 

The Coulomb interaction is of paramount importance in fractional QHE systems, where it gives rise to 
quasi-particles with fractional charge obeying anyonic statistics. It came as a surprise that e-e interaction 
plays a prominent role in integer QHE interferometers as well, even when their conductance is ~ e 2 /h so 
that the Coulomb blockade physics seems to be inessential. For instance, visibility in the MZIs and FPIs 
strongly depends on the source-drain voltage showing decaying oscillations, which have been termed "lobes" . 
The search for a resolution of this puzzle in the case of MZI has triggered a lot of attention 4^, 44 , 4|| 42 , 
4fj|, U HU . On the contrary, the extent of theoretical works on FPIs operating in the integer QHE regime 



is rather small [HJ H3, H3 . 

In this section we develop a capacitance model of the e-e interaction in a FPI and apply it to study 
the transport properties of the FPI in and out of equilibrium in the limit of weak backscattering. Our 
approach is inspired by the previous theoretical work [50] ■ Its essential idea is that a compressible Coulomb 
island can be formed in the center of the FPI between two constrictions (Fig. 1), which strongly affects 
Aharonov-Bohm oscillations. Starting from this model, we demonstrate that depending on the strength of 
the e-e interaction the FPI can fall into "Aharonov-Bohm" (AB) or "Coulomb-dominated" (CD) regimes 
observed in the experiments [2^, [29[. We also analyze the suppression of nonequilibrium AB oscillations 
with the increase of a source-drain voltage and find regions of both power-law and exponential decays, which 
explains experiments of Refs. 

The brief account on results of this section has been reported by two of us previously [lfj. Here we 
present technical details of our calculations and further elaborate on the qualitative picture of the interplay 
of interference and e-e interaction in the FPI which explains well the plethora of experimental data on the 
flux and gate periodicity of Aharonov-Bohm oscillations. 



4.1. Model 

We consider an electronic FPI of size L formed by a Hall bar with v edge channels and two constrictions 
(QPCs) that allow for electron backscattering between the innermost right/left moving edge channels with 
amplitudes T\(^ as shown in Fig. [6] (a). Right- and left- moving channels are connected to leads with different 
chemical potentials /i+ and respectively. In what follows, we take into account the backscattering in the 
lowest order, thus accounting for interference of maximally two different paths. For simplicity we assume 
the flight times along upper and lower arms (i.e. between two QPCs) to be the same, r + = t_ = r = Ljvp. 
We denote the magnetic flux threading the interferometer cell by <fi, i.e. an electron which encircles the cell 
once accumulates the phase 27r0/</>o, where (f>o = hc/\e\ is the flux quantum. 

The 2DEG in the QHE regime is divided into compressible and incompressible strips [20(. The filling 
factor in the n-th incompressible strip is integer. These strips are separated by much wider regions of 
compressible Hall liquid with a non-integer filling factor (compressible strips). The corresponding sketch of 
electron density profile p(y) in the FPI along y-axis is shown in Fig. [7] Let us denote by the boundaries 
between compressible and incompressible regions. Then a,k = Dk~ Vk 1S ^ ne width of the fc-th incompressible 
strip while bk — yZ-i ~ Vk 1S w idth of the fc-th compressible one. As it was shown in Ref. [20|, in the 
situation of gate- induced confinement of 2DEG in the QHE regime the widths bk 3> Xb, with Xb being 

1/2 

the magnetic length. At the same time dk scales as at ~ {bkXs) , so that in general the condition 
bk 3> dfc Xb is satisfied. In this picture compressible regions play the role of edge channels — the 
self-consistent electrostatic potential is constant through the compressible strips and can be controlled by 
connecting them to external leads. 
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We also assume that the filling fraction v$ in the center of the FPI exceeds v, giving rise to a com- 
pressible droplet (Coulomb island). The reason for that can be smooth (on a scale As) disorder potential 
fluctuations [21[. Let us denote by eNi the excess charge on the island (e < 0), with N$ being integer. On 
the scheme in Fig. [7] the boundary of the island is given by y~ . This value is quantized and changes abruptly 
when an electron tunnels between innermost compressible strips and the island through the incompressible 
strip. On the contrary, the boundaries of edge channels may change continuously because of a variation in 
external parameters, such as fi± and V g , or due to quantum fluctuations of electrostatic potentials on these 
compressible regions (see also discussion later). 




Figure 6: (a) Fabry-Perot interferometer with compressible island i. The innermost edge channels e are subject to backscattering 
at the QPCs; the remaining f T = u — 1 right(r)- and left(i)- moving channels are fully transmitted; their role is in screening 
of the interaction between the electrons of the channel e. Right- and left-moving channels are connected to reservoirs with 
different chemical potentials fi^ = eV v , rj = ±. Here, the gate g is depicted as a "plunger" gate, (b) Simple capacitance model, 
which uses geometry of the compressible regions. Right(r)- and left(i)-channels are each joined into one conductor with the 
widths b r and 6;, respectively. 




Figure 7: The structure of the edge states in the FPI. Cut through the center of an interferometer cell: (a) Landau level 
energies as the function of vertical coordinate y and (b) electron density p(y) (see also Fig.|BJ. Besides v completely filled LLs, 
the central region can sustain a partially filled Landau level whose occupation jVj changes by tunneling. Here we denoted by 
a k = y~l — J/jT th e width of the fc-th incompressible strip (shown by grey), which has the integer filling factor k. Compressible 
regions (white) represent the regions with non-integer filling fraction and have the width = y7_ 1 — y^ . The picture above 
corresponds to the case v = 2. 



Electrostatics of the FPI 

In the framework of the above model, we treat the e-e interaction in the FPI by using the constant 
interaction model with mutual capacitances C a p between four compressible regions — the interfering channel 
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(e); right- and left- moving fully transmitted channels (r, Z); the compressible island (i) — and the gate (g). 
These capacitances are denoted by C eg , C e i etc. We assume a large capacitance between counter- propagating 
innermost channels — thus they share the same electrostatic potential (p e — and also consider fx = v — 1 
right- and left-moving channels as joint conductors with potentials <p r (<pi)- Defining a capacitance matrix C 
with elements C aa — ^ 7 C Q7 -I- C ag , and C a p = —C a/ 3 (a ^ /3), where Greek indices span the set {e, r, I, i} 
and q a = —C ag V g is an offset charge on the a's conductor, the electrostatic energy reads 

£ = ^ Z (Q* - <7a) (C" 1 ) (Qp-qp). (88) 

Z * \ /a/3 

Total charge Q; = e (iVj + v<j)/(f)o) on the island is distributed on the highest partially filled Landau level 
(LL) and on v fully occupied underlying LLs (cf. Fig. [7]). Single electron tunneling is possible between 
interfering channels (e) and the island. We assume the rate of such tunneling process to be much smaller 
than all other energy scales in the problem, r e t h, hence Ni is quantized and is fixed for given external 
parameters (4>, V g , /J±). 

The mutual capacitances C a p can be estimated from geometrical considerations [22|. We regard the 
island as a disc of radius r and represent the compressible edge channels as concentric rings of the width 
b a and diameter L (here a = e,r,l) as depicted in Fig. [6](b). The edge channels are assumed to be thin, 
b a <C L. Therefore for estimation of capacitances we can neglect the difference between the radii of the 
island and those of edge channels, i.e. L ~ nr ~ 7ry^ : . A top gate, if present, is modeled by a plane situated 
at distance d from the 2DEG. Since the size of FPI cell is much larger than d, we treat Ci g , C eg and C r n\ g 
as a parallel-plate capacitors and find an estimate 



r 2 Lb e Lb r 

Cl9 ^ e 4d> Ge9 ^ e 2^d> Cr ^ e ^d> (89) 



where e = 12.6 is the dielectric constant for GaAs. The estimate for edge-to-edge (C er and C e i) and edge- 
to-island (C e i) capacitances can be found as a mutual capacitance of two conducting rings. In the limit 
k<(iwe obtain with logarithmic accuracy [23| 

^-^Inf^, C^lnf^. (90) 



7r \ a 



Here b r m — J^k 1 ^fc are ^ ne total widths of fully transmitted edge channels. Finding the mutual capacitance 
between a plunger gate and the island or the interfering edge channel is in general more difficult. Because 
of geometry, one can expect that Ci g and C eg in this case will be substantially smaller than the above 
estimate (|5§f for the case of a top gate. 

Let us now comment on the flux dependence of the electrostatic energy, Eq. (JSSJ). When the magnetic 
flux through the island is increased, 5(f) = irr 2 8B , the LLs are squeezed and the charge on the island (for a 
fixed boundary y~) varies as SQi = ev 5<fi/4>o- A similar effect of magnetic field on charges Q e ,r,i distributed 
on compressible circular strips is negligibly small because of the condition b <C r. Indeed, for a typical 
variation SB, such that &4>/4>q ~ 1, the corresponding modulation of these charges are 

^^fAW (9i) 



and we do not include them into Eq 



4-2. Results 

In this subsection we summarize our results and give their physical interpretation. The detailed derivation 
is presented in the next subsection. The qualitative behavior of the FPI crucially depends on the relative 
coupling strength of the interfering edge (e) to the fully transmitted channels (I, r), to the island, and 
to the gate. The essential parameters are the number of transmitted channels v* which screen the bare 
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e-e interaction in the interfering channel — as we demonstrate in the section [4.31 one has v* ~ 1 in the 
case of strong and v* ~ z/ in the case of weak inter-edge interaction — and the effective edge capacitance 
C e as defined below by Eq. There are also two characteristic energy scales in our problem: (i) 

charging energy Ec = e 2 /(7 e , or charge relaxation frequency wc = {v* /"k)Ec; and (ii) the Thouless energy 
e t h = t^ 1 = Vp/L. A relation between these two parameters depends essentially on the geometry of the 
experiment (most importantly, on the geometry of the gates). We will assume that the condition e t h <C wc 
is always satisfied, which simplifies a lot our subsequent calculations and enables us to get analytical results. 
This appears to be a proper assumption for most of available experiments. In particular, the value of the 
Thouless energy that can be deduced from the experiment of the Harvard group is e t h ~ 50 fiV [26] , whereas 
the charging energy is in the mV range |28j ] . 

4-2.1. Visibility, dephasing and the "lobe" structure 

In the limit of weak backscattering, rj <C 1, the differential conductance of the FPI, g = g mc + <7ab, is 
the sum of incoherent and coherent contributions. The incoherent contribution is 



z (V)=v-R u (eV)-R2.(eV), 



(92) 



where Rj* (e V) are the renormalized reflection coefficients (see Eq. (|96|) below) . The dependence of the AB 
conductance on external parameters — the gate voltage V g , the variation of the magnetic field AB and the 
bias V — factorizes into 

g AB (fx + ,fX-, ( f>,V g )=g(V)cos[cp AB (V g ,AB)}. (93) 
The AB phase <^ab will be discussed in the details shortly. The amplitude of the oscillations is 



g(V) = e- T / T *R 12 *(eV) 2 
with the nonequilibrium dephasing rate given by 



(i 



cos ( \eVr\ 



7T 

Av* 



\eV\ (RUeV) + R 2 *{eV)) - sin 2 

7T ZV 



(94) 



(95) 



In Eqs. (|92p and (|M|) we have introduced the renormalized reflection coefficients defined as 



R jm {eV)=R j 



UJ C 



eV 



l/u* 



r(2- i/v*) 



Rj*(£th) 



eth 



eV 



l/v* 



Rl2*{eV) = i?i2 \uj c t 



l/2v* 



eV 



i/2u* 2 1 / 2v *e 1 / v 
r(l - l/2v* 



= -Ri2*(eth) 



eth 



eV 



l/2i/* 



(96) 



Remarkably, in the last equation the amplitudes r\, r-i do not renormalize separately, rather the renormal- 
ization operates non- locally. A similar result was found for FPIs in the fractional QHE regime in Ref. (49l |. 
The relations are valid for bias in the range eth "C eV <C u>c- The above renormalization comes 
from virtual electron- hole excitations (being a precursor of weak Coulomb blockade 0, HBj]) and stops at 
eV ~ eth- O n the contrary, the dephasing rate r" 1 is caused by real e-h pairs excited by backscattered 
electrons and is proportional to the shot noise of the QPCs. There a simple linear dependence of the shot 
noise on voltage, which is valid in the absence of interaction, is modified because of the renormalization of 
reflection coefficients. 

A functional dependence on bias in the conductance amplitude (I94p stems from an oscillatory prefactor 
which has a characteristic scale 7reth- As a consequence, the amplitude or, equivalently, the visibility v(V) = 
\d(V)\/9inc{V) vanishes for certain equidistantly distributed values of bias. The resulting characteristic 
"lobe" structure of visibility is shown in Fig.[5]and is in agreement with experiments reported in Refs. (27l.[26j. 
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Figure 8: Total differential conductance g as a function of bias and magnetic flux (left panel), and visibility (right panel). 
Parameters are: v* = 2, uic T = 25, Hi*(e t h) = R2*(^th) = 0.2. 

4-2.2. Aharonov-Bohm oscillations 

In experiment one usually characterizes the FPI in terms of a pattern of its equilibrium conductance in 
the (B,V g ) - plane, which is governed by AB phase. We have identified four different regimes where the 
behavior of AB oscillations is qualitatively different (see Table [1} . In this table the parameter v* — the 
effective number of transmitted channels which screen the Coulomb interaction in the interfering channel — 
depends on the relative strength of the inter-edge e-e interaction. To distinguish between the limits of weak 
and strong e-e interaction we compare an inter-channel interaction energy ~ e 2 /C ea (here a — r,l) with a 
screened by the gate charging energy of the interfering edge itself given by ~ e 2 /C* ff . Here the edge-to- gate 
capacitance is effectively increased by the so-called "quantum capacitance" : 

C* ag = C ag + (y - l)re 2 /(M- (97) 

In the weak coupling limit one has C ea ^> C* ff . In this case the electrostatic potentials on all edge channels 
(r, I, and e) are approximately equal to each other and we set v* ~ v. In the opposite strong coupling limit 
we have C ea <C C* g . The potential <p e here fluctuates independently of potentials on other edge channels 
(r and I), thus screening of e-e interaction by the latter channels is not effective and one gets v* ~ 1. To 

Table 1: "Phase diagram" of the FPI, which discriminate between two Aharonov-Bohm (AB and AB*) and two Coulomb 
dominated (CD I and CD II) regimes. 





C ea ^> C* g 


C ea < C* ag 


V = V 
Cei = C e i ~\~ C r i + Cu 
Ceg = C e g + C r g + Clg 


V* = I 
Cei " C e i 
Ceg — Ceg 


C e i Cig 


AB 


AB* 


C e { 3> Cig 


CD II 


CD I 



make a distinction between the "Aharonov-Bohm" (AB) and "Coulomb-dominated" (CD) regimes we now 
define an effective edge-to-island capacitance C e i by the relation C e i = C e i + C r i + Cu in the weak coupling 
limit, i.e. at v* ~ v, and set it to be C e i = C e % in the opposite case of strong coupling. Then the FPI 
falls into AB or CD regimes depending on a ratio C e i / Ci g , as it is shown in the Table [TJ As one can see 
from Eqs. I|89p and ([TO)) , for the device with a top gate the capacitance Ci g scales like r 2 when the FPI 
size grows, while C e i increases only as rlnr. This suggests a simple rule of thumb: the AB regime occurs 
primary in large FPIs with a top gate (in experiment "large" means a cell area ^20/im 2 ). In this situation, 
i.e. at Cei I Ci g <§C 1, fluctuations of charge on the island are screened by the gate electrode and do not affect 
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the AB conductance. In the case of opposite ratio between the capacitances (C e i/Ci g ^> 1), as one will 
see shortly, the AB conductance becomes linked to the Coulomb blockade on the compressible island. This 
explains a terminology choice — "Coulomb-dominated" — for the above regime. 

For a device without the top gate a bare edge-to-gate capacitance C eg is due to only a plunger gate (see 
Fig. [5]). Such gate is used to control the size of the interference loop and because of geometry C eg typically 
is very small, so that one has C* g ~ [y — l)re 2 /(hw). In this case our first condition of weak versus strong 
inter-edge e-e coupling can be simplified. Defining the dimensionless coupling constant as 

_ (v-\)e 2 

a u = z 

envp 

and using the estimate (|90p for capacitances C er and C e ; one obtains the crossover value 

* 1 . fb„b r 
a ~ — In 
2tt 

which sets the boundary between the weak and strong coupling regimes. 

We name four regimes AB, AB*, CD I and CD II according to the Table Q] above which for the benefit 
of the reader lists the values of parameters v* , C e % and C eg . The capacitance C eg here is defined in analogy 
to C e i. In addition to these effective edge-to-island and edge-to-gate capacitances we now define full island 
and edge capacitances as 

C% = C e i + Ci g , C e — C e g + CeiCig/Ci. (98) 

Then the AB phase in Eq. reads 

^ab - - ^^(N„ + v4>/tM + M (2V g -V+-V-) (99) 

with integer iVj* minimizing the charging energy of the Coulomb island 

Ei = ^r(^» + Hl<t>o - C ig V g /\e\) 2 , (100) 

and we have defined u>c = v*Ecl^ and Eq — e 2 /C e - 

In Fig. [9] we show the conductance gAB(<t>,V g ) in the (B,Vg)-plane for three regimes: AB, CD-I and 
CD-II. The plots display significant differences. In particular, the lines of constant phase have a different 
slope in the AB and type-I CD regimes. The flux periodicity is also different in these two cases. The AB 
conductance in the case of type-II CD regime shows the "rhomb-like" pattern. The pattern of equilibrium 
conductance in the AB* case is the same as in the AB regime, provided one sets v* = 1. 

4.2.3. Discussion and comparison with experiment 

Let us now discuss the physics which underlies the rich phenomenology of our rather simple model. To 
appreciate the role of interaction we consider first the non-interacting case, where the interfering channel 
couples neither to the fully transmitted ones nor to the island. Consider electron contributing to the 
tunneling current which is, say, incident from the left source and leaks into the left drain. It may tunnel 
either at the left or right QPC. The latter path is longer than the first by 2L = 2vpT and encircles a magnetic 
flux 4>. Along this path the electron accumulates a dynamic ( "Fabry-Perot" ) phase 2er (e is the energy of 
the electron) and a magnetic AB phase 2-k4>/4>o- According to quantum mechanics, the current results from 
interference of both paths. Integration over all energies in the range /i_ < e < /i+ gives the back-scattering 
current 

- de 



h 



2tt 



ri _j_ e 27ri0/0 o +2ieT r2 



Iab, (101) 
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Figure 9: Aharonov-Bohm conductance: AB regime (left), u* = v = 2, stripes of constant conductance have a negative slope, 
flux and gate voltage periods are A<j> = <pQ and AV 9 = vEc/\e\\ type-I CD regime (middle), v* = 1, v = 3, stripes of constant 
conductance have now a positive slope, flux and gate voltage periods are A0 = 0o/(l — v) and AV g = Ec/\e\; type-I CD 
regime (right), u* = v = 2, Ci g /C e = 0.6, at fixed gate voltage V g conductance depends discontinuously on flux — phase 
jumps occurs at every A</> = <j>o/v. 

that we have split into incoherent and coherent contributions, 

L mc = -—V(R 1 + R 2 ), Iab = ^ ^cos[(M + +/i-)r-27r^ ]sm[eyr]. (102) 

Z7T Z7T T 

While the incoherent part of the current ii nc is expected already on the classical level, Jab stems from 
interference and is sensitive to magnetic flux. The dynamic phase accumulated by an electron depends on 
its "absolute" energy e, and hence the current Jab depends on both chemical potentials but not 

just on their difference eV = /i+ — Clearly, the sum (/Lt + + enters only into the phase shift of the AB 
pattern, but not in the amplitude. However, this independence of the amplitude of oscillations on the bias 
does not in general hold for the differential conductance $ab- Specifically, when the differential conductance 
is calculated in the framework of the model of non-interacting electrons, the amplitude of the corresponding 
AB oscillations g does depend on the manner in which bias is applied. 

Experimentally, the bias is applied asymmetrically: /i+ = eV , /i_ = 0. The expected conductance then 

is 

<?ab = G^OvIab = -2i?i2 cos[2eFr + 2^/0 o ], 

i.e. bias merely controls the phase shift of the AB oscillation pattern. The amplitude g — —2R\2 is 
independent of bias. This clearly contradicts to our results presented above as well as to experimental 
observations. 

The situation would change essentially if the bias were applied symmetrically: /i + = eV/2, ji- = —eV/2. 
Then the conductance would be 

9AB = -2i?i2 cos[eW] cos[27T0/M ( 103 ) 

Now, the amplitude would oscillate with bias on the scale 7reth, yielding a visibility with a "lobe" structure. 
This result is apparently much more similar to our findings (albeit without dephasing and renormalization 
of R12) as well to the experimental observations. On the basis of the similarity between Eq. (|103[) and the 
experimental observations it was conjectured in Ref. [2|| that the electron-electron interaction effectively 
symmetrizes the bias even if the latter is applied asymmetrically. 

To see how this works, assume that a charge within the interferometer cell produces a (for simplicity 
constant) self-consistent potential ipo. An electron which propagates in this potential during a time 2t 
accumulates the "electrostatic" AB phase — 2(poT. Hence, the dynamic phase would be 2(e— tpo)T and instead 
of the bare chemical potentials the relative potentials (/i± — (po) enter the result (|102p . Such a mean-field 
potential is indeed generated within our model. For instance, in the generic limit of large charging energy 
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Ec 3> £th our calculations in the subsection 14.31 yield ipo ~ +/x_)/2 in the case of AB regime. Therefore 
without a need of any fine tuning, the bias is effectively symmetrized, which explains the appearance of the 
"lobe" structure. 

The mean-field potential (po on the compressible strip corresponding to the interfering edge channel is 
in general the function of applied chemical potentials p± , the gate voltage V g and the magnetic flux <j). The 
most general expression found in the subsection 14.31 reads 



yo(M±,^ ff ,0) = - 



f 



+ UJ C T 



eV g + /i+ y~ c C T + E C ^{N„ + za^o) 



(104) 



Here Ni», as before, provides the minimum for the Coulomb energy Ei of the island given by Eq. (|100[) . If 
we introduce the electrochemical potential 



then the AB phase, given by Eq. (|5§]h is equivalently represented by relation 



VAB{v±,V g ,4>) = 2lT<f>/<j)o - 2Tip (n ± ,V & 



9' ' 



(105) 



(106) 



The first contribution here is the magnetic AB phase accumulated along a fixed reference loop with area 
Ao, i.e. </> = Aq5B, with SB <C B being a weak modulation of magnetic field on top of the high field B 
which drives the 2DEG into the QHE regime. Because of the condition l»« i (see Fig. [5] b) imposed in 
our model, and since a typical variation SB is such that <p changes on a scale of a few flux quanta only, one 
can use any boundary y k to define Aq. For example, one can set Aq = 7r(y+) 2 . Let us further show, that 
the second "electrostatic" contribution (— 2t4>q) to the phase (fAB can be interpreted in terms of a motion 
of edge states which leads to the variation of a relevant FPI area when the magnetic field B or gate voltage 
V g are varied. 

First, we note that in a stationary limit an imbalance of electron density per unit length on the interfering 
edge channel 8p e is related with the corresponding electrochemical potential f|105[) by a simple relation 
Sp e = —(Po/(2itvf), since (2ttvf)~ 1 is just but the ID thermodynamic density of states in our model. As it 
is always the case in QHE systems, this charge density can be translated into the variation h u of the boundary 
between the compressible and incompressible strips (Fig. [TO]), Sp e — n^h^, where = B/cj) Q = (27rA^) _1 is 
the electron concentration on one completely filled LL, and we have assumed the fluctuations of inner and 
outer boundaries of the compressible strip to be the same, Sy^ +1 = 6y~ = h v . Therefore "electrostatic" 
part of the AB phase reads 



2r(po = iTTL5p e = Ann L (Lh y ) = 2tt{B / (t> )5A. 



(107) 
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Figure 10: When the boundaries of the compressible strip (white) separating incompressible regions (gray) with filling factors 
v and v + 1 move towards the (y + l)-liquid (indicated by arrows), the latter shrinks while the i^-liquid's area grows: as the 
result the charge on the incompressible strip, defined by the boundaries y~, j and yt, decreases (in this figure h v < 0). 
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Here SA — 2Lh L , is the change in area enclosed by interfering edge state. To recapitulate the logic, we 
have thus related the self-consistent electrostatic potential ipo to a variation of the FPI area. With these 
arguments at hand one can now rewrite Eqs. (I104[) and (|106D in the equivalent form 

4>ab = 2n(A SB + BSA)/c/) + const, (108) 

where the variation of area reads 

^tB jB+ 1M + ^' (109) 

In the last equation we have introduced an integer which is a deviation in excess number of electrons 
on the Coulomb island with respect to the excess charge corresponding to some initially chosen reference 
gate voltage V g . We have also took into account the source-drain bias is small on a scale of the charging 
energy, \e\V < E c . 

Let us now discuss our theory of the FPI in relation to the recent experiments by considering separately 
each of the four regimes in the Table [TJ In the AB regime one has C e i/Ci <C 1, thus the coupling of the 
interfering edge to the Coulomb island is negligible and the area A does not change with B. The AB phase 
then simplifies to 

<^ AB = 2^/0 o + (110) 
v E c 

yielding the lines of constant phase with a negative slope (Fig. |H1 left) and a magnetic field period AB = 
(^o/Aq which is independent of v. The second term in the above equation describes a modulation in space 
of electron trajectory under variation of the gate voltage. If 5V g > then the interfering edge state moves 
outwards and thus encloses a larger flux as it is seen from Eq. f|109|) . The AB regime was observed in 
large devices (cell area ~18/xm 2 ) with a top gate (2(| [2|| [29|, where the condition C e i <C C g i is satisfied. 
In Ref. [28[ it has been also found that magnetic field period AB is independent of B but gate voltage 
periodicities (both top- and plunger- ones) scale as AV g cx l/B oc v. These observation are consistent 
with our Eq. (|110l) if the charging energy Ec = e 2 /C e is v - independent. It is so, provided the full edge 
capacitance C e ~ C eg + C e i stays approximately the same at each Hall plateau. 

In the CD regime one has C e i/Ci ~ 1 and the area of the interfering loop shrinks with the magnetic 
field. This is because interfering edge is now electrostatically coupled to the charge Qi — e(Ni + vcf)/(f>o) on 
the Coulomb island, which has explicit dependence on flux. The AB phase in this regime reads 

^ AB = 2tt(1 - v/v*)<p/J> - ^ < 277 Wa 



v v 
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In the type-II CD regime v* ~ v and at fixed V g the AB phase stays piecewise constant when the magnetic 
field is varied. The dependence of phase on B exclusively enters via 2Vj. Indeed, as it follows from Eq. (|109p 
the shrinkage of the interfering loop in this case, 6 A = —(Aq/B)SB, exactly compensates a change in 
magnetic phase tp = AoSB. When the FPI is brought close to a charge degeneracy point of the island 
by varying V g or B, electron tunneling becomes possible between the droplet and interfering channels (i.e. 
ANi = ±1) resulting in abrupt change of A. This creates a phase lapse (or jump) A(/>ab = ±2ir/v giving 
rise to the "rhomb-like" pattern shown in Fig. [S] (middle) at v > 2. 

In the type-I CD regime v* ~ 1 and a change in AB phase caused by area shrinkage when rising B 
overcompensates the magnetic AB phase, since now 5 A — ~v(Aq/ B)5B. Counterintuitively, the phase 
decreases for increasing the magnetic field. At the same time, whenever an electron tunnels into the island 
from the interfering edge channel (e) , the boundary of this edge state contracts so as to expel exactly one 
flux quantum from the AB loop. The phase lapse, being equal to A<j>AB = ( — 27r) in this tunneling process, 
is therefore invisible in the interference conductance. As the result, one has the diagonal stripe pattern 
with lines of constant phase having positive slope (Fig. right). The periods are AB = <f>o/(fxA) and 
AV g — e/(C eg + C g i), with fx = v — 1 being the number of fully transmitted edge channels (note, that at 
v = 1 the lines of constant phase are vertical). 
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In the limit of weak backscattering at QPCs, the Coulomb-dominated regime has been observed in 
Ref. 29]. In this work measurements were performed with the set of edge state configurations (including 
fractional fillings), classified by bulk filling fb and number fa of fully transmitted edges. We focus on the 
results obtained with a 4.4/^m 2 -device without a top (but with plunger-) gate. They indicate that interaction 
plays a major role (i.e. CD-I and CD-II regimes). For integer and Jt = /& — 1 the results coincide with 
our Fig ED (right), including the period in magnetic field which scales as AB oc l//x- The gate period in 
Ref. [29( was found to be weakly increasing with fr. We can explain this dependence if one assumes that the 
full edge capacitance, which equals to C e = C g i + C eg in the CD-I regime, decreases with fr, since a mutual 
coupling of the plunger-gate to the inner interfering edge channel (C eg ) becomes less efficient at high v. 
Quite "exotic" behavior was observed for more than one number of channels trapped in the interferometer 
cell. In case of fb = 4 and /t = 1 experimental findings resemble very much Fig. [5] (middle) corresponding 
to our CD-II regime. In such setting one fully transmitted and one partially reflected edge channel can 
be described by our model with v — 2 assuming that other two inner trapped channels play a role of a 
compressible island, where the excess charge is quantized. 

Non-equilibrium transport measurements in the FPIs in the AB regime have been performed in Refs. [26L 



27] . Their main findings can be summarized as follows: (i) the dependence of AB conductance versus B and 



the bias V factorizes into a product of two terms yielding a "checkerboard" pattern in the (SB, F)-plane 
(cf. Fig. [51 left); (ii) the scale of the "lobe" -structure is set by exh ~ Vd/L] (hi) a visibility decay with bias 
is stronger at higher magnetic fields. Results of our theory, Eqs. (|93p - (|96D . are in full accord with these 
observations. In particular, a suppression of visibility in our model at <C 1 is mainly due to a power-law 
decay with the negative exponent (— l/2z/*). Since in the AB regime v* — v oc \/B, this decay is stronger 
in case of a small number of edge channels, i.e. at higher B, in agreement with Ref. (26j . 

It is interesting to note, that our theory predicts a fourth regime (AB*, see Tabled]). It is characterized 
by the same equilibrium conductance pattern as the AB regime (Fig. [SI left), but in contrast to the latter, 
the power-law decay of the visibility oscillations corresponds to v* = 1, and thus is independent of B. Such 
a behavior of the FPI has not yet been observed in the experiment. 

Closing this section we have to mention that a crossover from the AB to the type-I CD regime (in our) 



terminology has been recently discussed in details in Ref. [5l|. We note that our capacitance model is very 
similar in spirit to the one used in that paper. However, the important difference is that our approach takes 
explicitly into account quantum corrections to classical geometrical capacitances, given by Eq. (|97[) . As the 
result we obtain the extra type-II CD regime which may arise because of screening of Coulomb interaction 
by the fully transmitted edge channels. 



4-3. Calculations 

We show here how to derive the above results using the formalism developed in Sect. [21 For simplicity we 
assume all edge channels to have the same length L and same velocity — vp (with > in right-moving 
channels and < in left-moving ones). Consequently, all flight times r = L/vp are the same. Further, 
the scatterer 1 has a coordinate x 1 and scatterer 2 has the coordinate x 2 = x 1 + L for each of the edges, see 
Fig.QH 
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We remind that two characteristic energy scales play an important role in our analysis, namely, Thouless 
energy e t h = t and charging energy Ec = e 2 /C e (or charge relaxation frequency u>c — ^-Ec), see 
Sec. sec:FPI-results. As has been discussed there, we will assume that the charging energy is much higher 
than the Thouless energy, and will consider the range of voltages intermediate between these two scales, 
uj c > |eV| > e th - 

4-3.1. Electrostatic Action 

Our network consists of v right-moving and v left-moving chiral channels which we label with index fx, 
see Fig. 1111 The innermost right-moving channel (p = +) is coupled to the innermost left-moving channel 
(p = — ) by two scatterers i = 1,2 with 2 x 2-scattering matrices s l . The remaining chiral channels (right- 
moving ones labeled \x — rl, . . . , r(v — 1), left-moving ones by fj, = 11, . . . , l(y — 1)) connect sources to drains 
without any possibility of tunneling. 

Interaction is taken into account by the electrostatic model described in the beginning. For simplicity 
we assumed that electrostatic coupling between all fully transmitted right-moving channels (p = rl, . . . , r[y— 
1)) is strong such that they share a common electrostatic potential V r . This enables us to merge them into 
one conductor (labeled a — r). We proceeded in the same way with the fully transmitted left- moving ones 
(p = 11, . . . , l{v — 1) are now merged into a — I) and the two innermost channels (// = +, — are merged into 
a = e). This reduces the number of charge degrees of freedom characterizing the edge channels down to 
three: 

Q e = dx^ + tp + + dx$-ip-, Q r = ^2 /dxVv K ?/v K , Qi = ^2 dz'tpiKipiK- 

As a fourth conducting element, we introduce a central compressible island with total charge Qi = e(Ni + 
v<P/4>q). While the second contribution, the charge in the v fully occupied LLs of the central region, is fixed 
by external parameters, the occupation N$ of the partially filled LL of the island is an (integer) degree of 
freedom to begin with. Since it is assumed to fluctuate via very slow tunneling, T <Si e t h, we will, however, 
treat it in a mean-field approximation. 

The fermionic action we start with reads as follows 

A [tp, $,Ni] = / dtdxipfi (id t +iv f _ l d x ) - {Q a - q a ) (c^ 1 ) {Qp-qp), 

where the first sum extends over the 2v chiral channels /i = rl, . . . , r(y — 1), 11, ... , l(y — 1), +, — and the 
second one over the 4 conductors a,(3 = e, r, I, i. The electrostatic part of ^4o[V'j ^i] 1S > 01 course, a direct 
consequence of ([55]) and we refer to the corresponding section for definitions of q a and C. Charges Q a are, 
of course, dynamic, i.e. time-dependent, quantities, but for the sake of readability we leave time-dependence 
implicit (as we did with time integration in the above electrostatic action). 

Our interest lies in interference effects which manifest themselves in tunneling corrections to current. 
Tunneling phases respond only to the Hubbard-Stratonovich field ip = eV e on the interfering edges (note 
that the electrostatic merging of channels fj, = +, — allows us to use just one field tp = ip + = if-). The 
short-term goal of the present section is to integrate out all other degrees of freedom. 

Potentials V c = (V e , V r , Vi) 1 : . First, we decouple the quadratic charge terms via a (multidimensional) 
Hubbard-Stratonovich transformation, thereby (re)introducing the potentials V a on the conductors. Since 
we do not need the potential Vi on the island, we single out the island degrees of freedom beforehand, writing 

G- X = (l cc P A, V ci =p c i(Qi-qi). (Ill) 

The index c refers to the 3 indices e, r, I, that means p cc , p C i = p\ c , pa are 3 x 3-, 3 x 1-, and 1 x 1-matrices. 
With that the action reads 

A in t[^,i>,Ni] = --(Qi - qifpu - -(Q c - qcfpcdQc - q c ) - V^(Q C - q c ) 
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and becomes upon Hubbard-Stratonovich decoupling: 

AntbP, & Vc, Ni\ - i(Q< - ft) 2 p M + ^(K - KOVccH^ - Ki) - V?(Qc - 9c) 

with y c = (v e ,v r ,ViY. 

Integrating out Q e , Q r , Qi: . Next, we integrate out the charges Q e , Q r , Q\. As explained in Sect. 12. ll 
and 12.21 charges and potentials are decoupled by a gauge transformation © which generates a tunneling 
term At (see below) and, according to the Dzyaloshinskii-Larkin theorem, quadratic and linear (in voltages) 
terms V 2 S a and Q a V a . The former is given by the polarization operators (fTO)) and amounts for screening, 
thus an "enhancement" of the capacitances (in fact, the capacitances become complex, "Keldysh"- and 
energy- dependent, but in the static limit the corrections are indeed positive). Then the retarded/advanced 
components of the "screening capacitances" read 

" r i _ P ±iuT 

S r J a {u) = -e 2 J2 Jd£ 1 dZ 2 Il r J°(w,Z 1 ,b) = ±i(v-l)e 2 — , a = r,l, 

~" * (112) 

1 - e ±lu)T 



Sll^u) = ~e 2 J2 I d6d6n^( W ,ei,6) = ±*2e ; 
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Charges injected from the reservoirs due to nonequilibrium boundary conditions (in excess of the equilibrium 
charge which is canceled by the positive background charge) are 

Q r = (y-l)e——, Qi = [v - l)e— — , Q e = e . (113) 

ZTT Z7T Z7T 

We collect them in the diagonal matrix S = diag (S e ,S r , Si) and the vector Q c = (Q e ,Q r , QiY- Subse- 
quent elimination of charge degrees of freedom transforms Aintfy, il>,V c , Ni] into 

MVcNi] = ~pti(Qi - ft) 2 + ^K'fc 1 + S)V C - V*(Q C - q c + P ~JV a ) 

with pu = p u - PicPccPd- (114) 

Integrating out V r , Vi: . The final, and somewhat cumbersome step is to integrate out the voltages V ri VJ. 
In order to do that we again split the degrees of freedom, writing 

s'-fe *)■ S =( S ' *)■ (115) 

where the index t refers to the 2 indices r, I, and c ee , c et = c* e , c u are corresponding 1 x 1-, 1 x 2-, and 
2 x 2-matrices. The action then reads 

A [V C , Ni] = - -piiiQi - ft) 2 + i(c ee + S e )V e 2 - V e {Q e -q e + C ee V ei + C et Vu) 

+ \ V t( C tt + St)V t - V t \Qt -q t + C te V et + C tt Vti - C te V e ). 



Performing the Gaussian integration over V r , V\ is a straightforward, albeit cumbersome calculation 
which in the end yields, 

Ao[V e ,Ni]=- i {Qi-q l ) 2 -{Qt-qt) t Pti{Qi-qi) + ^C*V 2 -V e {Q e -q e + x et (Qt-qt)+x e i(Qi-q i )) 

(116) 
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where we have introduced the effective capacitances and coupling strengths 

C' 1 EE p U -picPcc Pci + (PieCet + PitC t t){ctt + St)' 1 {c t ePei + CttPti), 
C* EE C ee + S e ~ C et {c t t + St^Cfe, C e EE C* 
X et = -Cet(c t t + St)' 1 , 

Xei = CeePei + CetPti ~ Cetffitt + S t )~ 1 (c t ePei + CttPti), 
Pti = (Ctt + St)~ 1 {c t ePei + CttPti)- 

Regimes ABI, GDI, CDII: . In the limits of very strong, i.e. 

C ae <C C e g + S e , 

and very weak, i.e. 

Q ^ [C ag + S a ) 2 

Ceg "i" Cg + Clg + S e + S r + Si 

coupling between fully transmitted edges a = r, I and interfering edge e the expressions simplify to 



c* e 

Xet 



Strong coupling 



c 



eg 



C e i + C, 
Se + C } 

(0 0) 

Cei I Ci 
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"/ — C ig /Ci 



Weak coupling 



Cei + C r i + Cli + Cig 
+ S r 



C e g + C r g + C; g + S e + 5/ + SV + Cig — Cf g / Ci 



(C e 



and using Eqs. (|lllj) - (IH5p . the action becomes 

1 



C r i + Cli) I Ci 

o, 



1 



= § j^tdt'ip{t)V-\t-t')y{t')- j^dt^(t)No(t)-^-(Qi- qi f 
with y r/a (w) = E c 

N 



\e\V g C e 



(117) 
(118) 



Ci 



■Qi/\e\ 



Here Ci, C e , C eg , C e i, and ^* are given in Sect. 14.2.21 fEa. ((98)) and table above). 
4-3.2. Tunneling Action 

To construct the tunneling action At[f] in lowest order we use Eq. (|2U|) which makes it necessary to 
identify the paths (ij;fj,v). Only the innermost chiral channels fj,, v = ± allow for tunneling between each 
other at scatterers i,j — 1,2 which gives 4 classes: (11; H — ), (22; H — ), (12; H — ), and (21; H — ). Classical 
phases are accumulated due to magnetic flux <\>: 



A^ 1 . = A<# 2 _ = 0, A<^ 2 _ = -A<j>f_ = -2ir(f)/(t) . 
At zero temperature the distribution functions read f±(t) = e~' iAt± * /q- (t) with Fermi distribution function 



fft-(t). Writing for short ee si 



+ • 



X = X+ - X-, £i. 



£ij3 (the right-hand side being the 3-dimensional 



Levi-Civita symbol), and Iljj ee we obtain the tunneling operators (|2"2l |2"3"]) 

n|(t) = -r^ e ±l * e -^[2^/0o+(M++M-)r] e -ie^ /o ^ (t + e . iT)/ § (t _ e .. r)) 

n£/ f (t) = [e(±t)u>(t) + e( T t)u<(t)} x ^. 
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(119) 
(120) 



Writing the tunneling phases $ = $^ the tunneling action reads 



-i ^2 dt i dt 2 (e~ l * e -** + (*V0 



n>(*12) nf(t 12 ) i l e i« + (^,ia) 



(121) 



with ti2 = ii — £2- According to Sect. I2.3l and Eq. @ the tunneling phases are related to the potential ip via 

&(x i ,t) = Q-(x i ,t)-Q+(x i ,t) = [dt'V + -{x t ;t,t')(p(t') (122) 



with V + -(x;t,t') = / dx'D -(x-x';t,t')- dx' D 0+ {x - x'; t, t') 



with bare particle-hole propagator Do/i (Eq. @). Note that x' is integrated over because potential (^(a;', t) ■ 
ip(t) in our model does not vary in space. 

Defining — ± for * = 1,2 retarded and advanced components of 2?_| read in energy representation 



V r /l(LU-X l )=±lE l 



CO 



(123) 



4-3.3. Current in Instanton Approximation 

Current is measured via the counting fields \ m the tunneling polarization operators (I119p . We use 
the adiabatic approximation where measuring time to is much larger than all intrinsic time scales of the 
system and transient effects due to switching of the counting fields are negligible. The tunneling correction 
to current is the derivative 



// = —i — <9 V In Z 



to 



I<-I>.) 

13 tj ) 



= E ft 

»,J=1,2 

with if = L Jdhdt 2 [vfYl e iA °&> N ^ +iA *M e -^ a (-\ti)nff(t 12 )e^^^ 



{Ni} 



(124) 
(125) 



The average J £VX){iv-} IS treated in real-time instanton approximation as outlined in Sect. 



e 

T 

e 



x =0, Ni=Ni, 



- dtidtee^^e-^^+t* ^ fijf (tia) 
to 



x =0 : Ni=Ni, 



(126) 



As always (. . .)o denotes averaging with respect to *4o[<P, ATj»] given in (|117l) . Because of the linear-in-tp con- 
tribution potential tp and hence tunneling phase $ have non- vanishing expectation values ipo = <po[Ni], $o = 
$o[A^i] which minimize «Ao[v> -^"i] f° r given iVj. At the saddle-point iVj* in turn minimizes ^4o[vo[A^], Ni]. 
For strong coupling u/qt 3> 1 the mean-field reads 



eV g + M+ ^ M u c t + E c ^r(N u + W/0o) 



1 + UJ C T 

$o(a; 1/2 ) = T^ =>■ -i<P (x k ) + i$ (x l ) = 2ie k iTip Q 



(127) 



with iVj, S Z minimizing the electrostatic energy Ei, (|100[) . 

Due to the presence of the source term i„4,/M = — i& a (x l , ti) + i&P(x J ,£2) the instanton phase 
<J>o + (|2"7|) . deviates from the mean- field by 



( xfc , £) = D^(t — ti,x k ,x i ) — D$(t — t 2 , x k ,x j ) 



(128) 
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The instanton action thus reads 

iAM = X)J>a /dt3dt 4 e- i * ( a *)^* ( a, )e- w *^»*' t »>+«»J( s, > t *>n2f(t3 -U) (129) 

kl -y5 J 

We will evaluate the time integrals in (|126p and (|129[) approximately. They will be dominated by 
the singularities of the instanton and the polarization operators. To identify and characterize them more 
precisely it is indispensable to compute the phase correlator D$ = —i (($ — < £ > o)( ( 3 ) — ^o))n- It wn l turn out 
that the singularities (branchcuts) of flfcfc(i) around t ~ and of !!&(£), k ^ I, around t ~ ±r dominate all 
integrals. 

4-3.4- Correlation Functions 

In this section we calculate the correlation function of the tunneling phases <I>(V) which according to 

(|122[) is D$ = — 2?_| VT^ Details of the calculation are not important for the rest of the paper and may 

be safely skipped. The final results for zero temperature and the strong coupling limit, loct S> 1, are 

D%(t = 0,x\x j ) = Dl /T (0,x\x 3 ) = D^-x\xi) = ( 7 + ln[w c r]) (130) 

and 



(t,x iX 1 ) — Aij 



In 



a ± i(t — t) 
a±it 



In 



for large times, \uct\ ^> 1. 

We start the computation by combining (111 8[) and (|123l) . 



Di{u,x uXj ) = -p;_( w; ^)y r H2?;_( W ;^') 



-iAi 



a±i(t + r) 

a at 



iuj c (l - e lUT ) 



(131) 



v* uj [u + iuj c (1 - e iuT )] 



with A^ = titj- In time representation the relevant correlation functions are the ^-components which, at 



zero temperature, read 

D%{t) = ± J 



A'jj 
2^ 



[DI{lj)-D%(u)]6(±w) 



[j*(t - r) - J*(t) + [j*(-t* - r)]* - [./*(-**)]*} 



with J^(t) = /dw [±0(±w)] 



io; c (l 



) 



cj [w + iuj c (1 - e iwr )] 



!)■ 



(132) 



(133) 



The integral defining J > (J < ) is perfectly convergent for all times with non-positive (non-negative) imaginary 
part, Imi < (Imi > 0), thus ensuring the analyticity of J< in this region. Apparently, we have J<(t)* = 

First, we perform the integration for Ret < 0. Under this assumption the contour of integration can be 
rotated into the upper half of the complex w-plane where the integrand is analytic (see Fig. [T2"]) . Defining 



dimensionless time and charging frequency, z 
imaginary axis, one obtains for y 3> 1 



-t/r, and y = u>ct, respectively, and integrating along the 



J*(t) = fd* (e— - 1) r USLlLJl „ [°° ds (e— - 1) -r^— 
Jo V ' s [s + y (1 - e- s )] J y > s(s + y) 

= - e y z r{0,yz)- 1 -lnyz = g(yz) (134) 
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Figure 12: Analytic structure of the integrand in (1 1 33 11 : it is analytic in the upper half of the aj-plane and possesses poles in 
the lower half. Contours of integration for J > and J < are indicated by arrows. 

with the incomplete Gamma function T(a,x) = J dse~ s s a ~ 1 , x € K, and the Euler-Mascheroni constant 
7- 

The asymptotic behavior of g is 



g{yz) =->(-i + i + ^yz)yz, 

-t - 7- Inyz , 

yz 



yz -> 0+, 
— > oo . 



(135) 



We now proceed with the case Re t > where the contour of integration can be rotated into the lower 
half of the complex plane cj-plane. In contrast to the previous case the integrand does possess poles in this 
region (see Fig. [T2|) . around which, therefore, the integral has to be taken additionally. Since both pole and 
imaginary axis contribution, and respectively, separately diverge for large lo we have to introduce an 
auxiliary ultraviolet cutoff, e^ aLJ , a = ar. Then, defining z = t/r, the imaginary axis contribution reads for 
y = io C T > 1 



> f°° 1 - e s 

Jo s [s -y{l- e s )\ 



1 e 



ds 



1 ±ish a±iz 
-e 1 = - In — - — . 



The poles are defined as roots of equation lu + iloc (l — e'" T ) = 0, lu ^ 0, and writing x — iujr, they are 
given by x n — y — W- n (ye v ), n g 22\{0}; where the product log function is defined by W n (x)e Wn ^ = x. We 
choose the numbering such that Imx„ + i > Imx„, Imxi > > Imx-i. As can be deduced already from 
the defining equation the roots satisfy Rex„ > 0. One may show that in two limiting cases one has 



x n —}2nin- 



V 



1 + y 
—}2win + In 



•21 2£ 
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Z7T 

w » h 

Z7T 



(136) 
(137) 



To proceed further, we note that \u> + iluq (l 
i"c (1 - e^) .. 



Res. 



LU [U! + iiOc (1 



1) 



1 + y — x n . Therefore the residues read 
1 



l + y -x n 



( e -zx n _ y 



Taking into account that for J> (J < ) only poles io n with positive (negative) real part contribute, n > 1 
(n < —1), we obtain for the pole contribution 



+ y- X n 



1 e 



±iax± n 
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This expression cannot be evaluated analytically further, but analytical approximations are possible by 
substituting the poles x n by their asymptotic behavior, Eqs. (|136l) . (|137l) . 

We convince ourselves that the short-time divergence, which forced us to introduce the ultraviolet cutoff 
a, is in fact merely an artifact of our method of calculation, and is cured by taking the sum of Jf+Jj. In 
other words, Jq- has to diverge logarithmically for z — > as well. Of course, any divergence originates from 
terms with large \n\ — > oo, such that for our present purpose we may safely use the approximation f|l 37[) 
which yields for z — » 



00 r / 



.2nnY 
V J 



1 , l-e- 2 ™ , ~a±iz 
— « — In . « In — - — . 

n 1 - e - 27r ( n± * z ) a ' 



which is exactly what we expected to find. Although the approximation is good enough to estimate the 

1 _ -2nd 

divergency, it is not reliable for obtaining finite offsets. Using In 5, = In - — |- — we can single out all 
cs-dependencies, 



4(t) 



Jf (t) = - In [±2mz] - J 



n=l 
00 r 

■E 

n=l 



1 



x± n =F 2-Kin — 1 — y 
(l + y-x± n )n 



-2iran 



The a-contribution is just a constant about which we will not care too much presently. For the moment we 
will fix it manually, by requiring a good agreement between Jq- (t) + J-f- (t) and the analytical continuation 
g{—ujc{t T iO)) of the result (|1 34[) obtained for Ret < 0. Fig. [T31 shows the corresponding plots for y = 5 
and y — 25. 




Figure 13: J > (real (blue) and imaginary (yellow) part), numerically evaluated and manually fixed, and analytical continuation 
~ g(—LOc(t — iO)) (real (red) and imaginary (green) part). 



A numerical study shows that the oscillating contributions decrease in width for large y (while their 
amplitude remains in the order of unity) and may be therefore neglected in the following. We approximate 
J< by smooth functions g<, required to be analytical for Imt < (Imt > 0) and 

(t) = - e-" ct r(0, -u) C t) - 7 - In [-u c t] , Re t < 0. 

Since the voltage is assumed to be low \eV\ -C ujc one needs correlation functions for long times |wc?*| S> 1 
only and we can use the asymptotic expression (|135|) for g. Therefore introducing a short-time cutoff a ~ ui^, 1 
and writing tip = t =p ia we use the following approximate relation in our subsequent analysis 

J^(t)w^(t) = -7-ln[-Wct T ]. 



which together with (fT32|) and J%(t = 0) = gives Eqs. (|T3T|) and ([T30]> . 
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Figure 14: Analytic structure of X in the complex t-plane: while the function is analytic in the lower half, it has poles or 
branchcuts in the upper half, residing at t* = 0, T, — T. Integrating n?-, (t) "around" these, i.e. along drawn contours, gives 



4-3.5. Renormalized Polarization Operators 

In the real-time instanton approximation, Sect. 12.41 virtual fluctuations around the instanton are taken 
into account by dressing the tunneling polarization operators, Eq. (pH]) . The phase factor is 



i[j3| oo (t,x»,x,)-D.(0)] =e ^( Wcr )^ 



a±it 




a ± it 


a ±i(t — r) 




a±i(i + r)_ 



and dressing of the bare polarization operators (|119[) yields (/o can be found in ([5]); x is put to 0) 



(2tt)2 



[a±it]^ 2 [a±i{t-T)] & [a±i{t + r)\ 
[a±it}~^ [a±i(t-T)]^ w ~ 1 [a±i(t + r)]- 



fc = /, 



The dressed polarization operators exhibit non-analytic behavior (poles or branchcuts) around t w 0, ±r. 



The double time-integrals (|126l) and (II 29[) can be approximately expressed in terms of the integrals P, 



kl — 

Jdt EK (t) . Before we demonstrate this statement in the next section, we will devote the remainder of the 
current section to the evaluation of Pj^ . 

We focus first on P^. To deal with both k — I and k ^ I simultaneously we generically consider the 
function 

X(t, eV) ee e- ieVt [a + it] v [a + i(t - r)] A [a + i(t + t)] A , 
0>?/>-2, > A > -1, 2X + n = -2. 

Apparently, I has branchcuts only in the upper half of the complex t-plane (Fig. IT4"]) . i.e. the integral 
JdtX(t, eV) vanishes whenever the integration contour can be closed in the lower half. Therefore, we assume 
the nontrivial case eV < 0. The real-time integrals J dtX(t, eV) = I_ +Xq +X + consist of three contributions 
which correspond to integrals along closed contours in the complex i-plane. With z = \eVr\ 3> 1 the contour 
integral around — r is 



2_ P3 - 2tt— — -e~ lz e^ x \eV\ Z - 2 - x . 
1 (-A) 

Similarly, one obtains the integral around +r, 1+ = X*_ . 

The situation is slightly less trivial for the integral around t ps 0, since it may be that 77 = — 1, i. e. we 
have a first order pole, or r\ > — 1, giving rise to a strong divergence. In the first case the integral gives 



Xq = 2n [a — ir] X [a + ir] X 



2ir 

T 
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In the second case we have to go around the singularity with care. We explicitly kept the distance i«a 
to the integration contour from the branchcut in the calculations. After weakening the degree of divergence 
by partial integration we may safely put 5 — > and obtain 



2A 



Zo w2 7 rr(-r ? )- 1 | e y|z 

Note that in this approximation Iq is continuous in n = —I. 

For the direct terms, k = I, we have i] = -p- — 2, A = — ^r, i. e. 77 + A = — 2, 2A = — hence for 
large voltages, z ^> 1, Iq is dominant. For the interference terms, k ^ Z, we set r\ = —-p-, A = — 1, i. e. 
77 + A = — — 1' 2A = — 2, hence the contributions I± dominate over Iq if and only if v* > |. 

As Jdtl(t,eV) splits into three contributions, so do p| = P,|(-t) + P^(0) + P^(+r). Note that 



n<(t)=n>(-t) 



-eV 



implies P<(f*) = P>(-f) 



-eV 



Summarizing, for z>l the dominant integrals of the dressed polarization operators are 

' ] 1 Rk*(eV), 





= 0) 




= kt) 




-0) 



v y 2tt 

0( Te y)_L R Wm (eV) ie lK [ |eV/T| -^( 1+ ^)]e- iefcip 

27TT 2 

0( T eV)-i- r k n \eVT\±- ^%]tf e^*** 

27TT T( — ) 



o+(M++j"-)t] 



(138) 
(139) 

(140) 



with k I, k — ±. In case of ^* > | the contribution (|139l) is dominant, while in the case 1/* < | it is 
(|140[) . We have used definitions of the renormalized reflection coefficients and assumed for simplicity 
r\fi to be real. 



4-3.6. Instanton Action and Current 

We have now everything in place to finalize the calculation of the instanton action (|129[) and the current 
(fT2"fJ|) . The instanton phases 8<5>Z{x k ,t) = Dl"(t - t 1 ,x k ,x i ) - D'f (t - t 2 ,x k ,x j ) and thus iA t [<p*\ are 
functions of the times t±, t 2 over which to integrate in (|126[) . A shift of integration variables £3/4 <-> t 3 / i + 1% 
in (|129[) immediately shows that the action i^4 t [y>»] is a function of the difference t = t\ — t 2 . Hence, the 
whole integrand of (I126p is purely a function of t = t\ — t 2 - Performing a change of integration variables 
(tiftz) i-> (t — t\ — t 2l T — {t\ + t 2 )/2), the integral over the center-of-mass time T is seemingly divergent. 
This simply amounts to infinite transferred charge Q = J *° dT I for a steady current I and an infinite 
measuring time to — > 00. Since our interest lies in the steady current (not on transient effects due to 
switching of the measuring device) we identify Jdtidt 2 — JdTdt = tg Jdt upon which the current becomes 



Iff = e dte iAt & 



(141) 



Given that iAt [</?*] is non-divergent, large contributions to this current stem from the singularities of 
fl^(i) which we identified in the previous section. The case v* = 1, i.e. v* < |, needs to be treated more 
carefully and will be considered towards the end of this section. Focussing for now on v* > 2, dominant 
contributions are then 



e 1 n 



e iA t [tp, 
piAt[tp* 



ti-t 2 f»0 



\tl-t 2 \~T 



I = J, 

r V0 (p«P( t * = +r) + P™f(t* = -r)) , i + J. 



(142) 



We evaluate the £3, ^-integrals in (11 29[) using a similar approximation scheme. Within the given con- 
straints t\ — t 2 fa for i = j and \ti — t 2 \ ps t for i 7^ j, the singularity of 



rifcfc(i3 — 44) 



(t 3 - k) 2 " 1 /"* 
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dominates over the ones of 



n kl (h-u)~- — — — .,1-1/2,. ' fc ^< 



and of the instantons 



e i&$>,{x k ,t' +t 2 ) _ e W^(t' -t 1 +t 2 ,x k ,x i )-iDq,(t' ,x k ,x j ) 



(t'-h+h) 2 \A hi /2u' / (t ,_ r)(t / +r) N^/2* 



(*' -tl+t 2 -T)(t> -ti+t 2 +T)J \ t' 2 

Hence, again transforming to relative and center-of-mass times, t = — £4, T = (£3 + t4)/2, the dominant 
contribution to the instanton action stems from the t ss O-singularity of Hkk (t) : 



k=l,2-y,6=T 

= - S A^fcfcW /dT (V J7 * (T) - l) with JT S (T) = -5$:(x fe ,T) + ^(x fc ,T) 



fc=l,2 7#<5 

The second equality follows from &>(*) + fi<(i) = ft^(i) + The integrals P^/ = Jdtn$(i) have 

been studied in the previous section. 



Using the definition (I128j) of the instanton and the relation (D^ — D^)(t,x k ,x l ) — (D^—D^)(t,x k ,x L ) 



Aki-p;0{—t)6{t + t) one obtains (independent of a and /3!) 

J*(T) - ±^ [A kl 9(h - T)6(T - *i +r) - A fei 0(i 2 - T)0(T - * 2 + r)] (143) 
v 

With constraints |ti — t 2 \ <^ t for j = j and \t\ — t 2 \ = t for i 7^ j (cf. Eq. (|142|l ). the instanton action thus 
reads 

iAAw } - I -l* 1 -* 2 !/^' , l ^ (144) 

t[¥M 1 -r/T v -ieVT€ ij (R lm (eV)-R^{eV))ir- 1 an^, i^j [ ' 

with dephasing rate 

V 1 = " 4shl2 2^ E (^(0) + P fc t(0)) (145) 

fc=l,2 

which due (| 138[) to becomes the expression (|95p in the limit \eVr\ 3> 1, see Fig. rj5] The purely imaginary 
contributions to (|144[) correspond to (perturbatively small) renormalization corrections to bias voltage and 
will be neglected further on. We will also neglect the instanton action for i = j. 

Combining these results with (|142D we obtain for the incoherent and interference corrections to current 
due to tunneling 

A/inc - hi + hi = ~ {Ri*{eV) + R 2 *{eV)) V, (146) 
Ztt 

Iab = -—Ru*(eV) signV e~ T sm{\eVr\ - n/Au*) cos tp AB (147) 

7TT 

with Aharonov-Bohm phase 

ipAB = 27r0/0o + (fjt + + ^_)t - 2rip Q (148) 
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Figure 15: Dephasing rate as a function of source-drain voltage shown for v* = 2, ioc T = 25 and -Ri*(eth) = R2*(^th) = 0.2. 
The solid line gives the numerical result for 11451 . The dashed line is the power-law asymptotic given by Eq. 1951 1. 



which in the limit ujct ^> 1 gives (jHU). For large bias, ujc 3> \ e V\ 3* £th, this yields the dimensionless 
conductances given in (|9"2|) and 

Concluding this section we turn to the case v* = 1. According to Eq. (|140j) the dominant contribution 
to current is 

e" 1 ^ w e ^dv.] e ? iet ' TVo P& , (t* = 0). (149) 

J tl-tjBSQ J 

The i j-contribution of Eq. (|142j) is also present, but subleading. The instanton action iAt [<p*] can be 
evaluated following the same line of reasoning as for v* > 2, yielding i_4.t[y>*] = — \t\ — t-^jr^ which can be 
neglected. The dominant contribution to current is thus 

I kb 1 = — — r l f 2 UJ C Te 7 cos^ab- 

TTT 

In the limit ojc t ^ 1 its bias dependence is negligible, in contrast to the contribution Jab, Eq. (|147[) . There- 
fore, while the latter is subleading in current for v* = 1, it yields the leading contribution to conductance: 
(Jab = dlpj^ 1 /dV + dlAB/dV w dlAB/dV, i.e. giving the previous result 



5. Summary and Outlook 

In this article, we have developed a general theoretical framework for investigation of electronic properties 
of a broad class of nonequilibrium nanostructures consisting of one-dimensional channels coupled by tunnel 
junctions and/or by impurity scattering. Our formalism is based on nonequilibrium version of functional 
bosonization. 

We have derived a nonequilibrium (Keldysh) action for this class of systems, Eq. (fl"2"|) . that has a form 
reminiscent of the theory of full counting statistics. In order to make the further analytical progress possible, 
we have developed a method based on a combination of a weak-tunneling approximation and an instanton 
(saddle-point) approach. 

We have supplemented a detailed exposition of the formalism by two important applications: (i) tunneling 
spectroscopy of a biased Luttinger liquid with an impurity, and (ii) nonequilibrium quantum Hall Fabry- 
Perot interferometry. One more application, quantum Hall Mach-Zehnder interferometry, has been presented 
recently by two of us with Schneider in a separate publication, Ref. [48( . 

The developed Keldysh-action formalism has allowed us to explore the rich interaction-induced physics 
of all the above problems, which includes, in particular, renormalization and dephasing far from equilibrium 
as well as an oscillatory voltage dependence of visibility of Aharonov-Bohm oscillations ("lobe structure"). 
The theoretical results are in good agreement with experiments on Fabry-Perot and Mach-Zehnder interfer- 
ometers. 

We close the paper by identifying future research perspectives; a work in some of these directions is 
currently underway. 



:-!<) 



(i) It is important to see under what conditions and how can one proceed in a controllable way 

(ii) A further important direction is the analysis of asymptotic behavior of relevant types of Fredholm 
determinants. Recently, such an analysis has been carried out for Toeplitz determinants arising in 
the problem of nonequilibrium Luttinger liquids and related models [53| . Required generalizations 
include, in particular, block Toeplitz determinants which arise naturally in the case of models with 
several channels coupled by tunneling. 

(iii) One of prospective applications of our formalism is related to edge states of quantum spin Hall topo- 
logical insulators. Of great interest is a generalization on setups based on fractional quantum Hall (or, 
more generally, fractional topological insulator) edge states. 
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Appendix A. Derivation of Keldysh Action 

The first two Appendices, |Appcndix A| and |Appcndix B| are devoted to the proof of (I12p and more 
specifically to the calculation of the Jacobian J[<p, ip,ip] corresponding to the gauge transformation ijy? — ► 

e ie ^ tp^ with (J2J. Since this is a linear transformation in ip, -0, its Jacobian is independent of tp, %p and can 
be therefore computed as 

J[<p] = e**tM = J V{i>,$)e iA ^^ = Tr {g U[<p + ]iU[<p-] } (A.l) 

with density operator go and the (many-particle) time evolution operators U{ip^]. The latter describes any 
single-particle dynamics the electrons undergo in the system and the leads (see Fig. IA.16|l . say scattering 
and propagation through time-dependent potentials <p a (t) . These potentials may have a non-trivial Keldysh 
structure, thus the superindex a — =F which refers to the forward/backward branch of the Keldysh contour, 
respectively. 




Figure A. 16: Sketch of the system (shaded blob) which is connected to reservoirs (rectangles) via source (drain) leads, depicted 
by incoming (outgoing) lines. The leads fi (u) enter (leave) the systems at contact positions x^ (xj?). 
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Appendix A.l. Integrating out Fermions with Klich's Formula 

There are several ways to evaluate the many-particle trace (|A.1[) . Here we employ an approach that 



In Appendix |Appcndix B we present an 
where an analogous action was derived for 



generalizes a derivation of the full counting statistics in Ref. 
alternative derivation which keeps closely to the spirit of Ref. 
the case of a single compact scatterer. 

A central mathematical statement proven in Ref. Q relates traces of certain many-particle operators 
with determinants of associated single-particle operators. We denote the (many-particle) Fock space repre- 
sentation of single-particle operators C by T(C) = cj (i\ C \j) Cj. Here, {\i)}i is some single-particle basis, 
and c, (cj) annihilates (creates) an electron in state Then, the following identity holds: 



det (l 



Tr I e 

To proceed we write the density operator in the form 

e r(F) 

Qo 



det (1 + e F ) 



(A.2) 



(A.3) 



where the single-particle operator F = ctiNi is a suitable linear combination of (single-particle!) "number 
operators" Ni — \i) (i\ in the reservoirs. E.g. in thermal equilibrium F — —fiHo with some appropriate 
Hamiltonian Hq — ^2 The many-particle time-evolution operator is canonically discretized as 



U[if a ] = Texp 



dtH[ip a (t)] 



N 



hm 

At-K> 

N-foo i=l 



-iAi r(H[ip a (ti)]) 



Hence, Eq. (|A.1|) is a trace over a (infinite) product of operator exponentials which, according to (|A.2[) . is 



J[<P] = 



Tr 



U[<p+]W[tp- 



det (t + e F U[(p+}^U[(p-}) 



det(l + e F ) 
= det (t + f (U[(p + ]^U[(p-} - l)) 



det(l 



(A.4) 



with the single-particle time-evolution operator J7[y Q ] (not to be confused with U[ip a ]) and the occupation 
number operator / = [l + e~ F ~\ 

Appendix A.2. Wave packet representation 

In a next step, we follow Landauer's original idea[55] and represent the time-evolution operators with 
respect to the wave packet bases, relating them to the single-particle scattering matrices. Using a more 
compact notation for the single-particle time-evolution operator 



U a {t u t ) = Texp 



i / dtH[<p a (t)} 



t < t u 



(hence C/[<^ Q ] = U a (oo, — oo)), Eq. (|A.4[) can be brought into the form 

JUp}= lim det{l + /[[/+(i + ,i_) t t/"(t + ,t_)-ll) 

f-{-— >-d=oo L L J J 



(A.5) 



We fix some time- independent reference Hamiltonian, say Ho = H[(p = 0], which contains the lead kine- 
matics as well as scattering, but no interaction or current counting. Then the incoming/outgoing scatter- 
ing states with respect to Ho form two natural bases of the single-particle Hilbert space, see Fig. IA.17I 
Each state is characterized by its energy e and the lead /i through which it enters/leaves the system: 



H M 



n/out 



:i/out 



The two bases are hence 



( M) 111 ) and (|e^) out ) e;AI . 
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Figure A. 17: (a) Spatial distribution of incoming scattering state \efi) ln (wavy lines). It extends in source lead fj, and arbitrary 
drain leads u, but in no other source leads, (b) Spatial distribution of outgoing scattering state |e!/) out (wavy lines). It extends 
in drain lead u and arbitrary source leads fi, but in no other drain leads. 

For the sake of argument we will assume the lead channels fj, to be one-dimensional (ID), semi-infinite 
and non-dispersive with constant velocity v^. We use the convention u M > 0, i.e. for source channels — oo < 
x < xf and for drain channels x® < x < oo, and choose the normalization such that ln /° ut (e'//|e/A ln ^ out — 
Sfj,fj,'6(e — e') is satisfied. 

The incoming state |e/i) m is a plane wave in source lead fi and spreads into drain leads v. It vanishes 
in all other source leads p! : 111 (x/i'|e/i) m = 8^ M -j= & ltx l v ^ (where \xfi) is the eigenstate of the position 

operator in channel fi). Analogous statements hold for the outgoing states. 

Let us now construct the incoming (outgoing) wave packet basis at reference time t- (t+). For that we 
define 

| f/i )-= [JLe^-*- x >^\e»), \t\) out = [^L e ^-^-^M\e\). 
J \/2tt J \/2tt 

Note that t is not a parameter which describes the time-evolution of a state "|^)" but labels the state |£/z) 
similar to e in |e^i). The two new bases are thus (|t/i) m )t ;M and (|i/i) out )t ;A1 . 

To shed light on the meaning of label t, we study the time-evolution of the newly constructed states with 
respect to the reference Hamiltonian Hq. For time t' one has 

in ( xu'le-W-*-) \tfi) in = V„ V^S(x -xS + v M (t - 0). 




Figure A. 18: Zoom into source lead /i: Sketch of spatial distribution of wave packet states \t-fi) ln , |ti/j) m , |t2A t ) m with 
t 2 > ti > t-. 

This is a wave packet in source channel \i which propagates towards contact x^, arriving there at time 
t, see Fig. IA.18I After entering the system it may split and spread in some complicated way. Similarly, the 
outgoing state |iA) out may be distributed in some complicated manner inside the system, however tuned 



42 



such that at time t it arrives at drain contact x® and continues propagation as a single wave packet in drain 
channel A. Summarizing, 

e -iH Q{ t-t.) m * = y/ ^\ x S) > e -iHo(t-t + ) | tA) out ^ ^ |^ (A 6) 

are wave packets residing at contacts x^, x® , and thus being independent of and t. 

Making the assumption that interaction, counting etc. is switched on and off adiabatically such that 
H{t') = Hq for t' $ [t-,t+], we now argue that the same simple relations hold when taking the full 
Hamiltonian Hit) into account, 

[/«(*,<_) M ln = V^|<>, u a (t,t+)\t\) m * = y fa\xZ). (Aj) 

These relations are a direct consequence of (|A.6[) and the fact that potential ip a is restricted spatially and 
temporally: For t > t' > t_ incoming wave packet U(t',t-) \t/i) is completely contained in the source lead 
where we assume ip a to be absent. For f < i' < t_ potential f a {t') is again ineffective since not switched on 
yet. Therefore, £/(*>*-) = e^^*"*-) |t/i) for all i. The reasoning is analogous for outgoing states. 

We are now able to give the operator /[/+(£+, t_)^[/~ (t + ,t_) in the wave-packet representation. For 
source channels /i, fi' the matrix elements read 

ia (tfi\fU+it + ,t^U-(t +) t^)\t'fi'f 

= E Jdt" U,»(tX)J2 /di''' in (^V'l^ + (i + ,i-)ni'''A) outout (i'''A|f/-(t + ,t-)KV) in . (A.8) 

^i" A 

Since the leads are populated by the reservoirs such that the occupation number of the incoming states is 
fixed, ln (efi\f |e'//) m = 8^^5(e — e')/ M (e), the distribution function in time domain simplifies to 

f^(t,t') = (t„| f |tV) = V / ^ e- fa C*-0 / M ( e ). 

The matrix elements of the time-evolution operator further reduce to 

out (a|t/(t + ,M |*V> in = out (t\\u(t + ,t)u(t,t')u(t',t_) |*» in = v^ov £/(*,*') |^) = ^(t.O 

(A.9) 

which defines the scattering matrix S — S [<p]. 

In the wave-packet representation Eq. (|A.5[) can be written 

J[(p] = e iA ' [ip] = dot [1 - / + 5 +t 5"/] 

where the determinant is to be taken with respect to source lead indices and arrival times. The log of this 
determinant appears in our general result stated in Eq. <\12\ . The retarded and advanced part of polarization 
operator which are present in Eq. (|12[) are not reproduced within the method of this section since they 
represent itself the quantum anomaly. Their structure does not depend on the actual noncquilibrium state 
of the system and can be deduced from the analysis of the fermion action in the absence of tunneling as it 
was discussed in section I^TTI 

Appendix A. 3. Construction of Scattering Matrix 

According to (jA.9|) the scattering matrix element S\^(t, t') is the transition amplitude for a peak residing 
in the (incoming) lead at time t' to a peak residing in the (outgoing) lead A at time t. One expects that it is 
the summed amplitude for all possible trajectories which connect /i with A. We will briefly demonstrate this 
assuming that the scatterer is a network of simple blocks which are connected to each other via "interface 
channels" which may be "outgoing" with respect to one block and "incoming" with respect to a neighboring 
one. The electronic state residing in the interface channel A is denoted by \x\). It corresponds to a wave 
packet which is leaving one of the blocks and about to enter another one. We further assume that the blocks 
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Figure A. 19: Sketch of the system: A block characterized by the dwell time r is connected to some part of the system (the 
shaded blob) via interface states \x v ). 



be simple enough such that each of them can be characterized by a unique dwell time r (possibly different 
for each block), i.e. a wave packet which enters the block at time t will definitely leave it at (exactly) t + r, 
through whatever channel: U(t + r, t) \xi) — Ufi(t) \x/) where |a^) is an incoming interface state and the 
sum extends over outgoing interface states \xf). This defines the functions u/i(t) and the scattering matrix 
elements 



where we have assigned a characteristic velocity v\ to each interface channel A. The full scattering matrix 
S\fi(t' ,t) can be constructed out of elements Sfi(t',t) with the use of following decomposition property. 
Consider the situation sketched at Fig. IA.19I Because of the decomposition property U(t',t) = U(t',t + 
r)U(t + T,t) of the time-evolution operator, the amplitude for the transition from a peak at t in channel \i 
to one at t' in A is 



Obviously, the inner transition amplitude can be decomposed further in the same manner, and the full 
scattering matrix element Sx^ turns out to be the sum of amplitudes , each of them corresponding to 
a possible path p connecting the incoming state \i with the outgoing one A. As p passes through a certain 
number of building blocks, is the product of the blocks' scattering matrix elements. For completeness 
we note that, since each trajectory will end in the outgoing leads, each decomposition will end with 



for A, A' £ out. 
Simple blocks 

Having convinced ourselves of the usefulness of definition (IA.10[) we turn to simple two examples of 
building blocks: wires with fluctuating potentials and point scatterers. Simple as they are, a broad class of 
devices, including quantum wire junctions and electronic interferometers, can be modeled as a network of 
these construction units, and in the following we will restrict ourselves to such systems. 

The corresponding scattering matrices are found by considering the time-evolution of wave packets 
tp(t,x) = (x\ U(t,to) \xi) which satisfy the Schrodinger equation 




(A.10) 



y/vxo^{x\\ U(t',t) |ag = /d<"V^V (*a| U{t',t") \x v ) s UIA (f,t) 



(A.11) 



V^a" (x x > I U(t', t) \x x ) = 5{t' - t)5 x >x 



(A.12) 



idtip(t) = H(t)ip(t), with initial condition ip(ta,x) = S(x — Xi). 



(A.13) 



We list the results here, giving all necessary definitions in the subsequent paragraphs: 
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Construction unit 




Scattering matrix 




Chiral wire x ; 


tp Xf 


M a (t',t) = e"^ (t ' 


>A(t',t) 


Point scatterer 


& 


Sf i (t',t) = s fi S(t' 


-*) 



Chiral wire with fluctuating potential. A possible block is a chiral non-dispersive wire where fermions prop- 
agate with constant velocity v in a fluctuating potential (p a (t,x). Similar to the leads wires are described 
by a single coordinate x, extending from xi to i/, such that the dwell time is r = — — — . It is taken into 
account by the "delay operator" A(t' , t) = 5{t' — t — r). The presence of the potential leads to accumulation 
of phase 



*/<(*) 



da;' tp(x' ,t — (x — x')/v). 



(A.14) 



The wire connects exactly one incoming to one outgoing channel and the scattering matrix has just one 
entry A4(t',t) in channel space. 

Point scatterer. Another possible construction unit is the point scatterer which connects one-dimensional 
incoming (index i) and outgoing (index /) channels such that scattering occurs instantaneously (dwell time 
t = 0+). The scatterer is characterized by the unitary time- independent scattering matrix s/;. If all 
channels A have a linear dispersion with constant velocity v\, then according to (|A.13I) the wave packet 
incident from channel i (at time to) is 



ip(t, x, A) = 5 X i S(xi + Vi(t - t ) - xj) + ^Sxf s fi J— S(x f + v f (t - t ) - x), t . 

f V Vi 

Since the state extends over several channels the wave function is a function of both channel index A and 
channel coordinate x, and the sum is to be taken over all outgoing channels /. 

Counting Fields. Up to now we have not addressed the issue of counting fields, claiming that they can be 
treated on the same footing as fluctuating potentials, a statement to be proven in this section. 

The number of electrons which flow through a certain point x in the time interval to < t < t is described 

quantum- mechanically by the operator N — f~dtl(t). The current operator / = vip^ (x)ip(x) becomes 

time-dependent in Heisenberg representation, I(t) = e lff ( t ~*-)ie~ lH (*~ t ~- ) (i_ is some reference time at 
which the initial state of the system is fixed), v is the fermion velocity in the considered channel. According 
to Levitov and Lesovik Q the correct generating functional of charge transfer through x is 



Z( x ) = (U- X (t, t ) f U x {t, to)) with U x (t, t ) = Texp 



i| dtl(t) 



From the properties U x (to,to) = 1 and idj:U x (t,to) = —\I(t)U x (i,to) we conclude that 



U x (t,i ) = e ^(*~-M e -*(ff-f/)(t-to) e -^(to-*-) 



is a possible alternative representation. Thus, 



z ix) = ( T e exp 



-i fdt'H x (t') 
Jc 



is the Keldysh partition sum with respect to the Hamiltonian H"(t) = H + v Jdx A^x, t)^ (x)ip(x) where 
time integration and ordering is to be understood along the Keldysh contour C and we defined the local 
vector potential A^(t, x) — a *ip-5(x — x) with the "time-dependent" counting field x(i') = X$(£~0^~ *o)- 
The corresponding scattering matrix reads s a (t',t) = 5(t' — t) e _J tx(t) anc i can b e incorporated in the 
total scattering matrix. 
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Appendix B. Alternative Derivation of the Action for ID Systems 

In this appendix we sketch an alternative derivation of the action A / [<p] which holds for one-dimensional 
(ID) systems. As shown in Fig. IB.20l thev may consist of several channels. Either direction of propagation 
is the same in all of them (in case of which the setup is referred to as "chiral") or there are two distinct 
possible directions: "right" (+) and "left" (-). The derivation generalizes that of Ref. Q and some of the 
arguments given already there will be not reiterated here. 




Figure B.20: Exemplary ID system with 3 right-moving (+) and 1 left-moving (— ) channels and 3 scattcrers. The system 
extends xl < x < xr. 



We use a single coordinate system xl < x < xr to describe all channels with left (right) contact position 
x l(r), i- e - velocities in right-(left-)moving channels are positive (negative): vf {—vf)- The fermionic action 

is 

AhM, V] = Jdx Jdt& (id t + iT 3 v F d x -E-V)^ 

where ^ = (tp") are vectors of Grassmannian fields with Keldysh and channel/direction indices a and /i, 
resp, t are Pauli matrices in direction space, Vf is Fermi velocity, £ is the self-energy correction due to 
the coupling to the reservoirs (see below), and V denotes the (temporally and spatially local) potential 
under the influence of which the electrons traverse the system. The latter may be the static scattering, 
Hubbard-Stratonovich or counting potential. 

To obtain the integrated action A, first its variation with respect to V is considered, 

iSA = - dx diTr r 



SV~(t, x)G T (x, t-0;x,t)- SV + (t, x)G T (x, t + 0;x, f)J = - Jdx Tr 5V(x)a 3 G(x, x) 

(B.l) 

where a are Pauli matrices in Keldysh space, and the first trace Tr r is taken over Keldysh and channel 
indices, while the "full" trace Tr is additionally taken over real time. The full fermionic Green's function G 
is needed at coinciding spatial coordinates, where it is discontinuous because of linear dispersion. The above 
time shift regularization, which takes into account that fermion fields in the action are normal-ordered, is 
equivalent to the identification 

G T ^ (x, t; x', t>) ^ J- (g T ^ (x; t,t')-S(t- t')) (B.2) 

with the quasiclassical Green's function 

g ap (x,t,t') =iv F (G a,3 (x + 0,x;t,t') + G at3 (x-0,x;t,t')) . 

Appendix B.l. Transfer matrices 

The Green's function G{x\,X2) is related to G(x' 11 x' 2 ) via the the single-particle transfer matrices 
M{x,x') for spatial evolution from x' to x: G a ^(xi,x 2 ) = M a (xi, x[)G alB (x[, x' 2 )M^(x' 2 , a; 2 ) t or, for 
quasiclassical Green's functions, g a ^(x) — Ai a (x, x')g a ^ (x')M^(x, x'Y . Defining U a (x;t,t') = ivp~ 1 5(t — 
t')Tz{id t ' — V a (x,t')) (or more clearly in energy representation U a (x)(e, e') = iwi?^ 1 r 3 (27r(5(e — e')e~ V a (x, e — 
e'))), the transfer matrices are given by 



M a {x,x' 



/s, _ / Xl exp dxi^ a (a;i)] , x >' x, 
d Bl exp[/*da:iW 01 (a;i)], x <' x, 
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where O x (O x ) orders subsequent operators with respect to their space coordinate x, smaller (larger) coor- 
dinates ordered to the right. Consequently, transfer matrices are diagonal in Keldysh space (with M. a being 
only related to V Q ) and satisfy 

5M a (x R ,x L ) = f " dxM a (x R ,x)[~iv F - 1 T 3 5V a (x)]M a (x,x L ). (B.3) 

J XL 

For short, we will also write for the total transfer matrix M. = M.(xr,xl)- Note that in chiral systems 
considered in Ref. Q scattering matrices S can be used. While coinciding with the transfer matrices in 
chiral systems, they differ in non-chiral ones, the two being related via 

S = r 3 (t+ + Mr-)' 1 (Mt+ + r_) r 3 (B.4) 



with the projectors t± = (tq ± t 3 )/2 in direction space. In chiral, say right-moving, systems T3 = t + = 1, 
r_ = and indeed S = M.. Generally, S is unitary, SS^ = 1, A4 is pseudo-unitary, A4(x, x')t 3 A4(x, x'y = 
T3. By defining g = o^gr^ spatial evolution amounts for the similarity transformation 

g(x) = Mix.x^gix^Mix^x'y 1 . 

The factor 03 ensures the normalization property g(x) 2 = 1. 

Using (IB.2|) . (|B.3|) and g(x) = M(x,XL)g(xL)M(x R ,XL)~ 1 M(x R7 x) the variation of the action (|B.1[) 
can be simplified to 

iSA = - [dxTr t 3 5V(x) — — (g(x) - cr 3 r 3 ) 
J 2tv F 

= ~\ JdxTr [-iv F ~ 1 r 3 5V(x)] M(x,x L )g(x L )M(x R ,x L )~ 1 M{x R ,x L ) + const. 

= - -Tr g(x L )M~ 1 SM + const. (B.5) 

where we absorb all contributions to the action which are independent of distribution functions in "const.". 
We will show later on that they vanish. 

Appendix B.2. Reservoir Green's Functions 

In a next step, g(xi,) is expressed in terms of the quasiclassical Green's function 

9L(R)[t ~ t] - Uw " w - 1 - 1 ( } 



of the left (right) reservoirs (with distribution functions fum which may have a non-trivial channel structure 
/i(jj)p). The Green's functions g~i — cr 3 giT 3 , i — L, i?, are related to each other via a similarity transformation 
as follows: First, introducing the Keldysh matrices 



one easily finds 



and hence 



.9, = UrVsTatfi with (7, = L> 4 i = 2 ^°) (B.<h 
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with = fi, = 1 — /j. Thus we haven proven 

9R = U~ 1 g L U with U = U£ 1 U R . (B.10) 

To relate the Green's function <?(xl) inside the system to its counterparts in the reservoirs one assumes 
that the dynamics inside the leads is governed by some relaxation process, say isotropization of momentum 
direction due to scattering off static white noise disorder. This is described by the self-energy contribution 



T,(t 1 ,x 1 ;t 2 ,x 2 ) = S(xi - x 2 ) 



Sfi(*i>*a)) xi,x 2 < x L 
2T re i/ \ 9R(h,t 2 ), x\,x 2 >xr 



Here T ro i denotes the relaxation time. The requirement that G{x + Ax, x) vanish for infinite distances Aa; 
yields the boundary conditions Q, 

(t+g L )(t-g(x-)) = 0, (l-g R )(t-g(x+))=0, (B.ll) 

again with g = a^gr^. Defining M. = UM. the second equation is equivalent to = (l — M^g^M) (1 + 
g{xh))- Combining it with the first equation gives 

= (21 + g L - M- x g L M) - (g L + M- 1 g L M) g(x L ) (B.f 2) 

and by inversion 

g{x L ) = 1 + 2(1 - g L ) (g L M + Mg^' 1 M (B.f 3) 

where we have made use of g\ — 1. 

To rewrite this expression we choose a specihc basis representation. Since g\ = 1 there exists one in 
which g~L = diag (1, —1). In the very same representation we write 



M 



Mu M12 
M 2 i M22 



Then we have gL-M+MgL — 2 diag {M\\, —M22), which is readily inverted, as well as 1+^l = 2 diag (1,0), 
l-g L =2 diag (0, 1), and (|B~T3f gives 

Dchning 

V=(t+g L )/2 + M(t-g L )/2 (B.f 5) 

one may show that 

g{x L ) =1-(1-, q L )V- x M (B.f 6) 

is equivalent to (|B.I4[) . We have thus expressed g(xi) entirely in terms of g~L(R) an d the transfer matrix 
M(x R ,x L ). 

Substituting this result into (|B.5|) and using Ai(x R , x£)~ 1 5M.{x R , Xl) = A4~ 1 SAi yields 



iSA = ^Tr (1 - g L )V- 1 5M + const. = TiV^SV + const. 

=S> iA = Tr Ln£> + const. = Tr Ln 



const. (B.f 7) 
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Appendix B. 3. Role of Drain Distribution Functions 

So far we have not been concerned with the channel structure explicitly and merely stated that there are 
2 directions of motion, right(+) and left ( — ), each of which is realized by a certain (not necessarily equal) 
number of channels (possibly even zero in chiral systems). In both left and right reservoirs to each channel 
\x was assigned a distribution function fe^ and For a right- (left-) moving channel (/l^) is the drain 
distribution function and one naturally wonders whether it should have an effect on the chiral fermions as 
long as the latter have not entered the drain reservoirs. We show here that this is not the case. 

To this end we introduce the block decomposition with respect to channel indices, e.g. 



M = 



M+- 
M- 



M+- 
M— 



fi = 



fi- 



i = L,R, 



(B.18) 



where e.g. M.^ , /;+, /j_ still may have channel structure (M.-\ )„„, (/»+)„) (/«-)„) however, with [i (v) 

extending exclusively over right- (left-) moving channels. 
Introducing 

= U„MU~ 1 -( Ql1 QA-(^M- + f<M+ -2f>M-f< + 2f<M+f>\ 

Q — UrMU l -l Qa Q22 )-[i M - + x M+ M -f<+M+f> ) (R19) 



and the projectors 



p + -i ( i +w ,)-fr °). p. -id- .*,)-(- I 



(B.20) 



the action (|B.17I) reads 

iA = In Dot (P + + QP_) + const. = InDct (P + + P-QP-) + const. 

where 



T+Q21T- T+Q22T4 



P-QP- 

Writing out the direction structure explicitly yields 



U-fcD* (fl-M^Mt- -^-^^^-^Uconst. 



iMX- M- ++ f< + + M + ++ f U 



This proves that the action depends only on and i.e. the source distribution functions. 



Appendix B.4- Tracing out Keldysh Structure 

We now return to the expression (|B.17I) . Since we already know that the result does not depend on 
fz,,- and Jr.+ we may make the choice /r,._ = and = In other words, we put gi, = 

gn = diag (g~L+, 9r-) = 9'm = osgmT^. As suggested by the subindex g ln contains the source (or incoming) 
distribution functions. It can be parametrized analogously to (IF3.9|) . 



(B.21) 



1 {2f > — 2 f < 

^3.9in = U^asUin with U in = -j= fj^ _ J t 



with / being the matrix of source distribution distribution functions. Since U is not needed anymore and 
hence M = M. , we recall the definition of V given by Eq. (IB.15|) and of scattering matrix S (|B.4[) and obtain 



V 



1 + 



M- 



,1 



OS 9T3 



(t + +Mt-)t 3 U^ 



1 



0-3 



UinSUr 1 



1 



(B.22) 
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With 



Q = u in s ur 1 = 



Q- 



Q- + 
Q ++ 



f>s- 



-ls- 



f<S+ -2f>S-f< + 2f<S-f> 
f- ±5+ S~f< + S+f> 



the action (|B.17|) is 

iA = In Det V + const. = In Det ( 



const. = TrLn [l - / + 5 +t 5"/] + const. 



(B.23) 



Up to anomalous terms, representing the R- and A-parts of the polarization operator, and contributions 
which are independent of distribution functions /, this proves Eq. (|12p 

Appendix B.5. Constant contributions 

What are the "const. "-contributions to the action we keep ignoring? All that we know so far about them 
is their independence of distribution function /. According to (|B.23j) they can be recovered by substituting 
/ = in the full action (where "const." is not neglected). So let us put / = in the rest of this section 
and recalculate the full action. Eq. (|B.13|) can be evaluated explicitly now, yielding g T (xL) — 9q{xl) and 

9 T (xl) = 9^{ x lV witn 

/ x ( 1 0\ ^ „ f 0^ 



1 



1 



where A vll (x;t' ,t) is defined as the amplitude for a /i-wave packet at position x and time t to end up as a 
f-wave packet at time t' and the same position x (cf. Sect. jAppcndix A.3[ ). The above relation between 
g T t T and holds for all positions x with g^{x) = Ai a (x, xi,)go (xL)M. a (x, xlY- To calculate the latter 
we define (for given x) the "scattering matrix" 

s = t 3 (t+ + M(x, x l )t-)~~ (M(x, x l )t + + T_) T 3 
of the region between xl and x which takes into account paths extending only within this region. Using the 

M(x,x L ) M(x R ,x) 




Figure B.21: Spatial evolution of g from xl to x with _M(a;,2;£); s is the corresponding "scattering matrix"; the shaded blob 
represents the rest of the system. In contrast to s, paths contributing to A ^(x) may extend throughout the entire system. 



recursion relations (see Fig. IB.21[) 



A ++ (x) 



-A-+{x), A+_(x) 



- s__^_ + (x)s ++ , 
+ s+-^-+(x)s + _ 



A (x) = v4_ + (x)s + _ 



9o(x) = 1 



(B.24) 



one obtains 

.(x) A + ^(xY 
.(x) A__(x) y 

For a smooth potential V(x) all paths contributing to A ufJ: (x) have a non- vanishing flight time such that 
A y ^(x]t,t) — 0. At point scatterers, which are a mere idealization, the Green's function g(x) is not well- 
defined and we stay sufficiently far away from them. Concluding, for equal times the Green's function 
(|B.2[) vanishes. Substituting this result into Eq. (|B.f [) yields SA\f=o = 0, hence const. = ^4|/=o is indeed a 
physically irrelevant constant. 
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Appendix C. Weak Tunneling Regularization 



In general, the scattering matrix S a = <S[y a ] depend on time in a complicated manner which makes the 
exact evaluation of the functional determinant (TT21 unfeasible. In this section we present an approximation 
scheme which applies for systems with weak tunneling. A tunneling action At is derived by subtracting a 
clean action An from the full one A. By construction Ab is small and may be therefore expanded. The 
crucial step, the separation of actions, can be performed at zero temperature where reservoirs may have 
differing chemical potentials. As it will turn out, interaction effectively leads to a dressing of the point 
scatterers by phases $. 



Appendix C.l. Clean Limit 

The clean limit was discussed in Sect. 12.11 The action An is obtained from ([12")) by replacing S by the 
clean scattering matrix S" = diag (e l 1 A±, 



e I1? «Ajv) = & %i> A with the delay operators A M and the 



phases accumulated along the complete paths — : 
The main step in the quest of the tunneling action 



,.D 



At = A — An = iTr Ln 



[l-f + s+^s-f) (1 - / + S+t£-/) r 



is the inversion of the second bracket, a problem which is already dealt with in Ref. |54| . Since S 1 * is diagonal 
in channel space, the operator can be inverted for each channel separately. So we consider 

(i - f + sps-f)^ = t-U + a;VK~-Oa^ 

where A M (t',i) = 5(t' — t — t^) is an appropriate time delay operator. Note that it commutes with the 
stationary distribution function f^. For zero temperature we follow Ref.[H3], according to which the inversion 
problem can be reformulated in terms of a certain Riemann-Hilbert problem, which is solved by defining 
the phases and by 



with t± = t±iO. They satisfy V 



t± - if 



idt' [B(t - T 5(t - f )] [#-(t>) - 0+(tO] , (C.l) 
0+, hence d~ - i9 A = d+ - tf v = 0, and /„e ! ^/„ = f 



ffJ. e ^ fn — /p 



Thus, one obtains 



e-<{t-U)+e-<U 



■A 



and, taking now the full channel structure into account, 



At = — iTrLn 



with the "regularized" scattering matrix 



(C.2) 



Defining in energy representation N(u>) = —9(~uj) as the zero temperature limit of the Bose distribution 
function, we note that the function B introduced in (jC.ll) is the "generalized" distribution function B(ui) = 
1 + 2N(uj) ubiquitous in bosonic Keldysh formalism. This fact allows us to endow the phases with a Keldysh 
structure. We address this issue in the next section before turning to the computation of the regularized 
scattering matrix in the subsequent one. 
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Appendix C.2. Phases and Keldysh Structure 

Before we turn to the computation of the regularized scattering matrix let us consider the Keldysh 
structure of the phases d. To this end we consider the clean path x^ — > x^ containing the point x. An 
electron propagating from the source reservoir to x (arrival time t) accumulates the phase 



0ZM,x) = -v- 1 / dx'^ a (x',t-(x-x')/v^ = -D r ^ a Jx,t); 



(C.3) 



likewise when traveling from x (departure time t) to the drain reservoir x® it accumulates 



(C.4) 



dx V a (x, t + (x'~ x)/v^) = D^(x, t) 

where we introduced the retarded/advanced bare electron-hole pair propagators 

Dl'*{t- x', x) = ±6{±t)5(x' -(x + V)) 
along the the clean path. Quite obviously, the phase accumulated along the complete clean path satisfies 
W = ^l- m ^^-{x^-x)/v tl )^dl out {x 1 t-{x°-x)lv,) = - {D r 0fi - D^) ^{x,t-{x^-x)/v^). (C.5) 

The propagators directly make reference to propagation velocity and thus encode spectral (kine- 

matic) properties of the electron-hole pairs. As usual they lack information about the system's state which 
one expects to be contained in the missing Keldysh component Dq^ of the propagator. In this sense the 
phases $ are related to the potential (p via the linear operator Dq which does not have the full Keldysh 
structure. 

To overcome this "deficiency" we recall the definition made in the previous section 



1(B + !),?-- I 



(C.6) 



Combining now (|C5[) and (|C6[) and using the stationarity of B(t,t') = B^{t — t') we obtain 

+ (x° - x)/ Vfi ) = -l(B + l) (D r 0fl D a J tp-(x,t) + \{B - 1) (D^ D^) <p+(x,t) (C.7) 
With that the phase that an electron, traveling along a piece of wire between x\ and #2, accumulates is 

^21 (*) = " ^.out(^) + ^.outC* 1 ,*" ~ ^)Av) 

^ej(.T 2 ,<)-ej(x 1 ,t-( a; 2 -x 1 )/« Al ) 



with the phases 



-^.«wt(^*) + ^/*(* + (^-«)/«p) 

-\ [(B + t)D^ -(BT + \ [(B - t)D^ -(B T t)D^] <p+(x,t) 



(C.8) 



We thus have managed to rewrite the scattering matrix of chiral wires A4 a (x2,xi) = e 1 ® 21 A(x2, x\) — 
e 10 "^ 2 ) A(x2, iije' 9 "' 11 ' in terms of phases with full Keldysh structure. Formally we started with expres- 
sions which did not contain any distribution functions whatsoever and somewhat artificially included them 
by redefining phases. The usefulness of such construction will become apparent in the next section. 
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Appendix C.3. Construction of "Regularized" Scattering Matrix 

We have previously shown that the full scattering matrix S can be constructed out of simpler units. We 
prove in this subsection that the same statement holds for the regularized scattering matrix S. 

Since the full scattering matrix elements are amplitude sums S"^ — J2 p Au$ a over all paths p connecting 
x^ with , definition (1C.2[) directly implies that — Y) Ay^ is the sum of the regularized amplitudes 

= e-^A^A-'e-^A^ = e-^^A^e m ti<) (C.9) 

over the same paths p with phases 0^(„) along the clean paths x^ — > x® and xf, — > x„ . Note that 
ft^out( x v >*) — an d ^out( x ^ t) = &u(i + ( x u ~ x h)/ v h)j an d thus definition (|C8[) implies 0J(x^,£) = 

We now consider a generic path p which starts at incoming lead channel fi, winds through an alternating 
sequence of wires and point scatterers (numbered 1, 2, . . . , N with positions x%,X2, ■ ■ • , a; at), eventually ends 
in outgoing channel v' . The corresponding full amplitude is 

< )T = M*(x?, x N )s» x ■ ■ ■ s%M* {x^x^M^x 1 , xl) 
and simply becomes 

A%* = A v {x^ x N )e-^(* N h» x ■ ■ ■ e^^A K {x\ x^^s^^A^ 1 ^) 



upon regularization. This implies that regularization of the amplitudes amounts to regularization of the 
scattering matrices of the building blocks: The effect of fluctuating potentials (and counting fields) is 
removed from all wires and incorporated in the tunneling phases 

^{xi,t) = &l{x\t)-Q*{xi,t) 

which dress the point scattering amplitudes s{, . This result is summarized in the table below: 



Construction unit 
Chiral wire 
Point scatterer 



Regularized scattering matrix 
A(t',t) 

S^(t',t) = s^ jW*) 8{t'-t) 



Appendix D. Second Order Expansion 

We prove the second order expansion ([2~T]) in tunneling strength ( "tun" ) starting from (|20p . Obviously, 
Avfj, = Y2x(^XvV^xu * s ^ ne sum over all paths winding forward and along the forward Keldysh branch 
(i.e. with potentials (p~) from source [i to some drain A and then backwards along the backward branch 
(with potentials to source v (for short: "source fj, — » drain A — » source u"). The backward part p + is 
obtained by time-reversal of a physical, forward path p + (source v 
has amplitude A^ x = A x v v ^ (note that hermitian conjugation, f , reverses the order of partial amplitudes in 



drain A) with amplitude A Xfi 



and 



i(p) 



a(p + )U(p-) 



the product). Denoting with p~ the forward part of p = p~ ®p + , the total amplitude is Ay^ 

In the chosen order of accuracy we only take paths with a total of 2 or less tunneling events into account 
The expanded tunneling action is At — A\ + Ai with 



Ax =iTr ^(1 - A)^ = i ^ Tr 



A 2 ^-TrJ2 A ^ A ^ = lT, E Tr^PfyA^U 



(D.l) 
(D.2) 
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The sum in (ID.1[) extends over paths p±; "source \i — > drain v — > source fi" with 2 or less tunneling 
events (obviously the number has to be even), while the sum in (ID.2[) extends over paths p 2 — p 2 ® p'i'- 

" source \x — > drain n -^ source v — > drain A -i> source \i " , 

P2 Pa 

or equivalently, " source v — > drain A — ^ source [i — > drain k — > source v " 

p'2 p'2 

where the equivalence is ensured by the cyclic invariance of the trace. Since fj, 7^ v subpaths p' 2 and p' 2 ' each 
contain at least one tunneling event, i.e. in our approximation exactly one. Depending on whether tunneling 
occurs on the forward or backward Keldysh branch k, A = fi or v. It is quite obvious that no matter how 
often a path evolves in time forward and backward (once for Ai, twice for ^2), as long as it starts and ends 
in the same reservoir and contains exactly 2 tunneling events, it exactly involves 2 different channels /i 7^ v. 
We will use this fact for a systematic classification of all paths. 

To warm up we consider first the simplest (and least interesting) paths, which contribute to A\ and are 
of the form "source /1 — > drain /1 — > source /i" without any tunneling taking place. Time delay operators are 
canceled exactly (since forward and backward paths coincide geometrically) and no phases are accumulated 
(since after regularization they are carried only by the tunneling amplitudes). The total amplitude is thus 
a product of forward scattering amplitudes s^„ of scatterers j along the clean path — > : 

j j \ ^P- j J vj^y, 

All other relevant paths contain exactly 2 tunneling events. Not surprisingly there is a whole plethora of 
them and we are well advised to proceed systematically. To this end and according to the observation made 
before we classify these paths with respect to the pair (fi, v) of different channels fi 7^ v involved and the 
two scatterers i and j (possibly i = j) at which tunneling /1 —¥ v and v — ¥ [i respectively occurs. Note that 
in this classification classes (ij; \iv) and (ji, v\i) are identical. What classes (ij; [iv) are possible, of course, 
depends on the topology of the considered network (since e.g. not all scatterers are even connected to a 
given channel /1). But, as we will show, once a class is fixed, its contribution to At is essentially independent 
of topology! 
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Figure D.22: The 5 tope-logically distinct configurations for class (ij;fii/). All channels except for fi,u and all scatterers except 
for i,j ("distorted white circles") are dropped. Scatterers i ^ j are different in A-D which allows for 4 different orderings along 
channels fj, and i/. In E tunneling occurs twice at the same scatterer i = j. 



Fig. ID.221 shows the 5 topologically distinct configurations of channels /i, v and scatterers j, j for given 
class (ij;fiv). Arrows at the scatterers indicate direction of tunneling. All scatterers except for i, j are 
dropped. They are involved only in forward scattering events and thus lead to corrections C(tun 3 ). In this 
way propagation from source x^ to scatterer i just leads to an amplitude Ai.^.in which accounts for the 
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finite flight time r* in . Analogously we define Aj.^.u, and the same for and further 

y~ = A -1 s i_ A- • — A^ 1 p-^fi 31 ') A- ■ h { 
v+ = A" 1 (V+^A- ■ - A -1 A. ■ ¥ 

X+ = A" 1 (V'+^A- ■ - A" 1 e^O^A . 

For paths of given type p\ or pi fixing the Keldysh branches on which the 2 tunneling events n — > v (at i), 
v — > [i (at j) take place characterizes the corresponding paths uniquely provided they exist. Whether they 
exist or not depends on topology. Remarkably, they always do if the 2 tunneling events are required to occur 
on different Keldysh branches. Fig. ID.231 shows the paths for i : \i v occurring on the forward, j : v — > /i 




Figure D.23: Paths for topology A with i : fi — > v on the forward, j : v — > fi on the backward branch. The forward parts are 
represented by solid thick lines; the backward parts by dashed lines. The first path is of type pi, the second one of type p2- 

on the backward branch for configuration A. Similar paths can be drawn for all other configurations. As an 
example, we consider the first path in Fig. ID.231 which is of type p\. It starts at source n, tunnels at i on 
the forward branch, arrives at drain is, evolves backwards, tunnels at j on the backward branch and returns 
to source /i. Its amplitude is A^^ = X^X^. The contribution of both paths to At is 

where tunneling polarization operators Hij-^u are defined in (|22j) and we have made use of — — P vil + 
C(tun 2 ). Similarly •A^j. flv with Hfj v is obtained by considering paths with tunneling at i on the backward, 
at j on the forward branch. As mentioned previously the very same paths can be drawn for all topologies 
and yield the same action. 

The story is less neat if both tunneling events occur on the same Keldysh branch: While p2-type contri- 
butions to Afj.^ and Ajj.^v are still independent from topology, pi-type contributions are not as universal. 
Taking all paths carefully into account results in the table below 





H 7 - ;(J y(£)/ [ S Ufj, S tti\ 




A 
B 
C 
D 
E 


f~Z(t)fc(-t) + f<W>(-t) + f<(t)f<(-t) 

f<(t)f>(-t) 

f > (t)f<(-t) 
-f?WH-t) + hihit) + f<(-t))s(t) 


ft(t)f<(-t) + f<(t)f>(-t) + /<(*)/<(-*) 

f<(t)f>(-t) 

-fZW>(-t) + !(/<(*) + /<(-*))*(*) 



with f£(t) = fj?( t + t* in - in ), = ff(t— ii in + tI in ). We have incorporated the tunneling-free 

contribution (|D.3|) into the case E which amounts to the 5(i)-terms. 

Quite miraculously, using the symmetry relation f^~/ v (t) — ~f^/ v (t) f° r t ^ 0, one can show that f|23[) 

TIT 

holds, i.e. Ily'^ can be represented in a form which is independent of topology. We exemplify the proof, 



e -i<VO') n < e i* M „ 
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which is a straightforward calculation, on configuration A. In this situation, in > r*, ir , and r* iT , > rf, 
and thus />(*) = -/<(*) for i > and />(-*) = -/<(-*) for i < 0, hence, 



fi.m 



-f<(t)f<(-t) = o(t)f>(t)f<(-t) + e(-t)f<(t)f>(-t), 
ft(t)f<(-t) + f<(t)f>(-t) + f<(t)f<(-t) = 9(-t)j>(t)f<(-t) + e(t)f<(t)f>(-t), 

proving our statement. In the case E it is necessary to drop a constant (albeit infinite) and thus physically 
irrelevant contribution to the action to obtain (|23l) or (|24p. (See similar discussion after Eq. (|24p ) 

Appendix E. Saddle-Point Approximation 



We prove here the formulas of Sect. 12.41 To keep things readable we resort on a rather symbolic notation, 
in which the contributions to the exponent in (|26p read 



Ao[<p] = -^pV 1( P~ l P9o, A t [ 



-le 



-**MlW*( 2 ), Aj[p] = -Jp 



with the tunneling phases $ = —T>p being linear functionals of ip. The saddle-points of Aq + At + A,j and 
Aq + A j are denoted t^** and (/?*. We show first that up to the chosen accuracy C(tun 2 ) the former is not 
needed and calculation of the latter is sufficient, using that tp* — p** = 0(tun 2 ). A Gaussian expansion of 
the full action around p** reads 

{Aq + A t + Aj) [p] PS {Aq + A t + Aj) [<?**] + ^{p - 99**) S 2 {Aq + A t + Aj) [p„] {p - ¥?**). (E.l) 

Note that the first order term vanishes due to p** being the full saddle-point. We now successively replace 9?** 
by p* in the above expression. Expansion around <p** and using 6 2 (Ao +Aj)[p**] — V~ x — 5 2 (Aq + Aj)[(p#] 
yields 

(Aq +A t + Aj) [p^ = (Aq + A t + Aj) [<£>**] + C(tun 4 ), 
6 2 {A +A t + Aj)[p4 = 5 2 {A + A t + Aj)[<p„] + 0(tun 4 ). 

Writing p — p^ — (p — p*) + {p* — 95**) and using again p>* — p^ = 0(tun 2 ) Eq. (|E. 1[) thus becomes 
{Aq +A t + Aj)[p] PS {Aq +A t + Aj)[p*} + -{fp - p*) 5 2 {Aq +A t + Aj) [p*] {ip - <£*) 

+ (p* - 93**) 5 2 {A Q + At + Aj)[p*} {ip-<p**). (E.2) 

Only the last term, which is linear in p — p*, contains To proceed, we define h = {p* — 99**) 5 2 {Aq + 
At + Aj)[p*] = 0(tun 2 ) and perform the functional integration ([2l)| in Gaussian approximation (|E.2|) : 



Vve i{Ao+A t +Aj)[ v ] „ e i(Ao+A t +Aj)[<p,] (Det 6 2 {A +A t +Aj)[<p*\) 



-1/2 



x exp[-- h S 2 {A Q + A t + Aj)[p>*] h] 
2 *- .. ' 



ze i(A +Aj)[ v ,] exp 



=C(tun 4 ) 

iAM - ^TrLn (1 + V5 2 A t [pA) 



(E.3) 



Hence, p* dropped out completely. The first exponential satisfies e l (Ao+Aj)[ip t ] _ ^ e *Aj[<p]^ where (. . .)o 
denotes the average with respect to the free action .AqM- This is shown using (tp)a — p = V qq and 
({ip ~ <p){p — p))a — iV. Since -4q[v?] is Gaussian we have the simple relation 



e iAj ^ Q = {e- Uv ) =exp 



-iJ<P~ \ ([J{tp-<f)} 2 )) 



g — i J<£> — -| JVJ 
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On the other hand it is i(Ao + Aj)[tp*] — — |(^o + J)V(qo + J) = — § JVJ — iJip — ^qoVqo- he. up to the 
last term — ^QqVqq, which is canceled by normalization, we have 

e i(,Ao+Aj){<p*] _ e -ij<p-±JVJ _ ^ e iAj[<p]^ 

which proves ((29)) . 

To deal with the second exponential in (|E.3|) we expand the logarithm to leading order in V (again, a 
Gaussian expansion), 

iA t [<p*\ - rTrLn (l + V6 2 A t [v>*}) « iA t [<p*\ - -Tr V5 2 A t [<p*\ = iA t [<p*} + i (6A t [<p*](<p - <p)} 

z Z s v / 

=0 

+^ ((v - ^)<5 2 A[^*]((^ - tp)) Q « i (Ab - + V*])o • 
Defining the phases $ = —"Dip, $„ = —T)(p* it is 

(A^ - ^ + ^]> = -i (e-WMW+MDl n 12 e <[*(a)-*(a)+*.(3)]^ = _ ie -».(i) n 12 e ^ 2 \ 

i.e. Eq. (f30|) . where IT12 = e^ jD ' I '^ 1,2 ' )_ ' D ' , '^ 0,0 ' ) - > ni2 are the renormalized tunneling polarization operators with 
the phase-phase correlator 

2) = -t (($(1) - $(1))($(2) - $(2))) Q = (WD) (1, 2). 
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